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BACKGROUND OF THE INVENTION 

1 . Field of the Invention 

[0003] The present invention relates to a method for analyzing spectra including contributions from 
scattering so-called inverse scattering analysis using a renormalized form of the Lippmann- 
Schwinger equations, and to a system implemented on a computer and attached to an analytical 
instrument where a spectrum of interest is received by the instrument and analyzed in the computer 
using the renormalized form of the Lippman-Schwinger equations of this invention. 
[0004] More particularly, the present invention relates to a method for analyzing spectra including 
contributions from scattering so-called inverse scattering analysis using a renormalized form of the 
Lippmann-Schwinger equations, where the renormaUzed equation permits absolute and uniform 
convergence of the equation regardless of the strength of interaction in the system from which the 
spectrum was obtained, and to a system implemented on a computer and attached to an analytical 
instrument where a spectrum of interest is received by the instrument and analyzed in the computer 
using the renormalized form of the Lippman-Schwinger equations of this invention. 

2. Description of the Related Art 
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[0005] Many spectral characterization include inverse scattering components resulting from internal 
reflections of an incident waveform. These inverse scattering components can give information on 
both near field and far field properties of the object being analyzed. However, traditional 
application of the Lippmann-Sch winger equations to analyzed spectra including inverse scattering 
components a re 1 ess t he s atisfactory b ecause t he L ippmaim-Schwinger e quations o ften d o n ot 
converge or give oscillatory solutions that must be truncated to product approximate and sometimes 
misleading analyses. 

[0006] Thus, there is a need in the art for an improved mathematical theory for analyzing inverse 
scattering components that always permits solutions because the equations absolutely and imiformly 



converge. 



SUMMARY OF THE INVENTION 



[0007] The present invention provides a method for analyzing inverse scattering components of a 
spectrum of an object of interest, where the method utilizes equations that are absolutely and 
uniformly convergence and amenable efficient iterative computational determination, with leading 
terms allowing for fast tentative identification of the object from which the spectrum is obtained. 
[0008] The present invention also provides a method for analyzing inverse scattering components 
of a spectrum of an object of interest, where the method utilizes equations that are absolutely and 
uniformly convergence and amenable efficient iterative computational determination, with leading 
terms allowing for fast tentative identification of the object from which the spectrum is obtained, 
where the method involves obtaining a reflectance and/or transmission spectra of an object of 
interest using an incident waveform (electromagnetic or sonic). The spectra is then analyzed using 
the inverse scattering equations of this invention implemented on or in a processing unit (digital or 
analog) to derive a potential fimction representing the object. Generally, an adequate potential 
fimction can be derived from the first few leading terms of the iterative solution of the equations, 
where few means the first four terms, preferably the first three terms and particularly the first two 
terms. 

[0009] The present invention also provides an analytical instrument including an excitation source 
for producing an incident waveform, a detector for receiving either a transmission spectrum or a 
reflectance spectrum or both a transmission spectrum and a reflectance spectrum of an object or 
volume of interest, and a processing unit for analyzing the spectra, where the processing unit 
includes software encoding the inverse scattering method of this invention. 
[00010] The present invention also provides a sonic analytical instrument including a sonic 
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excitation source for producing an incident sonic waveform, a detector for receiving either a sonic 
transmission spectrum or a sonic reflectance spectrum or both a sonic transmission spectrum and 
a sonic reflectance spectrum of an object or volume of interest, and a processing unit for analyzing 
the sonic spectra, where the p rocessing unit includes software encoding the inverse scattering 
method of this invention. 

[0011] The present invention also provides an electromagnetic analytical instrument including an 
electromagnetic excitation source for producing an incident electromagnetic waveform, a detector 
for receiving either an electromagnetic transmission spectrum or an electromagnetic reflectance 
spectrum or both an electromagnetic transmission spectrum and an electromagnetic reflectance 
spectrum of an object or volume of interest, and a processing unit for analyzing the electromagnetic 
spectra, where the processing unit includes software encoding the inverse scattering method of this 
invention. Of course, the analj^ical instrument can include both sonic and electromagnetic 
excitation sources and detectors. 

We present a new inverse scattering series for quantum elastic scattering in three spherical 
dimensions. The new series, which converges absolutely, independent of the strength of the 
scattering interaction, results from a renormalization transformation of the Lippmann-Schwinger 
Fredholm integral equation to a Volterra form. A new feature of the formulation is that it does not 
require determination of phase shifts, and it can be applied even to integral cross section 
measurements. The approach is illustrated by application to a simple example problem. 

DETAILED DESCRIPTION OF THE INVENTION 
[0012] The inventor has found that the inverse scattering problem encountered in many field of 
spectroscopy can be solved in a more rigorous manner by applying a simple but effective 
renormalization condition on the traditional Lippmann-Schwinger equation. The simple and 
effective remomalization of the Lippmann-Schwinger equation results in a set of equations that are 
absolutely and uniformally convergent for all systems regardless of the strength of the interactions 
between the probing waveform and the obj ect being probed obviating the convergence problems that 
plague traditional Lippmann-Schwinger analyses. In fact, the inventors have found that calculation 
of only a limited number of terms provides an accurate enough representation of the object being 
probed, either near or far field, for rapid identification. The renormalized equations find application 
in all electromagnetic and sonic spectrometry application and is especially well-suited for analysis 
of data from very low frequency electromagnetic imaging used to tract undersea objects such as 
submarines. 
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[0013] In this application, we introduce an approach to the inverse scattering series that completely 

solves the problem of convergence[7 A] . This is achieved by renormalizing the Lippmann-Schwinger 

equation from a Fredholm to a Volterra structure. It was proved that the resulting inverse Bom series 

converges absolutely and uniformly independent of the strength of the interaction. However, the 

issue of how best to deal with the need for half-o-shell scattering information remains outstanding. 

We b ase o ur s trategy f or s olving t he i nverse p roblem one xploiting the superior c onvergence 

properties of the Volterra-based inverse series in combination with the introduction of a 

parameterization of the interaction that allows for a greatly simplified determination of the 

interaction parameters. This enables us to take account of half-o-shell (near field) effects. 

INVERSE SCATTERING THEORY: STRATEGIES BASED ON 
THE VOLTERRA INVERSE SERIES FOR ACOUSTIC SCATTERING 

I. INTRODUCTION 

[00141 It is also important to note, and we emphasize the fact, that use of a Volterra-based inversion 
requires either a) measurement of both the reflection and transmission amplitudes in order to achieve 
the simplest form of inversion b) the solution of more complicated, nonlinear algebraic equations 
if only the reflection amplitude is measured. In Section n of this apphcation, we present a simple 
analysis that shows clearly the need to deal with half-o-shell or near-field effects in order to treat 
the problem. In Section EI, we present a general analysis of the moments of the Fredholm and 
Volterra Bom approximations compared to the moments of the true interaction. This shows that the 
Voherra-based Bom expansion yields one higher moment before the half-off-shell effects come into 
play. Then in Section IV, we consider inversion of the scattering produced by any interaction that 
can be expressed as a sum of Dirac delta fiinctions (a model interaction that is shown in more detail 
in Appendix B to be of practical utility). This interaction also nicely illustrates the convergence 
properties of the Volterra series. For simplicity, we restrict ourselves to a scalar scattering wave (z.e., 
acoustic scattering), but the analysis is valid also for more complicated electromagnetic scattering. 
In Section V, we discuss the data required for implementation of the Volterra-based inversion. 
Finally, in Section VI, we give our conclusions. 
IL Analysis of the Role of Half-Off-Shell Transition Amplitudes 

[0015] The Bom-type inverse scattering series is most simply obtained from the abstract Lippmann- 
Schwinger equation 

T.^V^VGIJ, (1) 
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where V is the (local) interaction, 7^ is the transition operator, and is the free (unperturbed) 
causal Green's operator, 



(2) 



modified for the acoustic (and electromagnetic) case by the multiplicative factor, 1^, This factor 
arises because, unlike for the quantum scattering case, the scattering interaction is of the form ifc^F. 
This extra factor, ^, causes important changes in the scattering behavior compared to the quantal 
case as is elaborated in the Appendix A. Here 1^ is essentially the square of the spatial wavenumber 
parameter, describes the unperturbed wave propagation. Also, we explicitly indicate that the 
abstract operators 7^ and depend on A: as a parameter. The scattering amplitude for the process 



k -^k' then is given (using Dirac notation for compactness) 



(3) 



[0016] Taking advantage of the local character of V, we then have that 



(A:'M^>=-5-fJ^e'(*-*'>T(z) 



(4) 



[0017] It is important to note that the matrix element </c'|F|A) of the true local inter-action only 
depends on the difference (Ac - k'). Conversely, \i{k'\V\l^ is only a function of (A - k") then Kis 
local. If Vi^i) is real, then we also have that 

{k\V^ = {k\V\k'r (5) 
It can then be established that {k\V\^ and {k^T^ cannot simultaneously depend only on the 
difference {k - k) for if that were the case, then Equation (1) could be written as 

{k]T,\k) = {k'\V\k) + ldk"V{k"-k'){k"\G;,\k")T,ik- k") 



= {k^V\k)+ jdyVik- k'-y){k"\G^,\k")f,(y) 
= {k'\V\k)+ jdyVik-k'-y) 



(6) 



{2k - y)y + is 



-L(y) 
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[0018] The last line, which results from substitution of the Green 's function of Equation (2) into the 
previous line, clearly shows that the integral, and hence the LHS of this equation, is not solely a 
function ofk -k\ which is a contradiction. For inversion, we write Equation (1) as an equation for 

{k'\y\k) 



Using the resolution of the identity 



1= j^/r|rXr 



this equation becomes 



/ > , \ / \ f {k'\V\k"){k''\T,\k) 



(7) 



(8) 



(9) 



which is an exact equation. We first observe that if we take |A:'| = \k\ so that {k'\Ti^\k) corresponds to 
an on-shell matrix element, then </c'|F|A") still involves both on-shell and half-on-shell T matrix 
elements. 

[0019] As we have seen, although V \s taken to be a local operator, it cannot be true in general that 
is local (or, more precisely, (A:'|7;|^> is not solely a function of A - A:'). The approach of Moses, 
Razavy, andProsser [3A-5A] involves additional expansions of{k^Vlc) and (A^ ^TJi) in a power series 
of the on-shell (reflection) amplitude. Thus, in addition to the issue of convergence of the Bom 
expansion of Equation (1), there is also the question of convergence of these expansions solely in 
terms of the far field reflection amplitude. We note that this analysis also applies to the Volterra 

inverse integral equation. One can, in general, write the Volterra Green's operator, Gn. , as 



Ok 



plus a solution of the homogeneous free Green's function [7A, 8A] equation Such homogeneous 
solutions can always be written as a sum of separable, totally on-shell operators having the form 



Ok ^ll\(l>nk){Xnk 



(10) 



n 



where 



ik'-H,)\<f>^)={e-H,)\x„,)=Q 



(11) 



Then, 
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(12) 



and. 



= vG,j, + voj, = v[u o,rJ+ vg,j, 



(13) 



[0020] We define by 



(14) 



where, 



(15) 



[0021] Again one can show that 7^ is, in general, non-local, and a parallel analysis to that for 

holds. The major distinction between using Equation (3) to generate an inverse series for 
and using similar matrix elements of Equation (14), expressed 



is that the iterative solution of Equation (14) for F is guaranteed to converge absolutely and 
uniformly no matter how strong Fis [7 A, 8A]. There remains the question of whether the general 
matrix elements of V and can be expanded in convergent series of far- field amplitudes. Thus, 

whether one bases an inversion on or , both require the equivalent of half-on-shell 

information for their implementation. We next examine the behavior of these two inversion 
alternatives with regard to their relations to moments of the true interaction, because this shows an 
important distinction in how far-field quantities affect these two. 
III. ANALYSIS OF MOMENTS OF THE INTERACTION 
[0022] The moments of the interaction are defined as. 



as 



(16) 




(17) 



[0023] It is clear from equations (7) and (16) that V(n) is exactly given by 




(18) 
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V[n] = f J dzz" f J dke-^"^ [ ( 



-km 



(19) 



or 



(20) 



[0024] It is generally assumed that (-A|r4|A)and ^- kT^ k) are obtained experimentally by a far- 



field measurement. The terms involving ( - k 



VG* T 



k ) and the analogous matrix element of Tu 



are the source of far-field effects in the above equations. Again, defining, 



(21) 



and. 



V,{z) = [yke-^'^l^- k\T,\k) 



(22) 



we have, 



V{n-\= V,[n]- rdzz" f^J/te""*'/- k 

v— GO w— 00 \ 



vg;,t, k 



(23) 



and, 



V[n] =V^[n]- [yzz" £J i^yfee"''*^ (- k\VG,,f,\k 



(24) 



[0025] We shall now prove that, in general, V^[n] is exact (i.e., V[n] = V^[n] through « = 1 
while Vi[n] is only correct for « = 0. To do this, we substitute into these equations the coordinate 



representation matrix elements of Gq^ and Gq^ , which are given respectively by [8A] 



- ik 



(25) 



and 
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(z|z') = k sm(k[z'- z^H(z'- z) 



(26) 



Here H{z ' - z) is the Heaviside function,. 



H{z) = 



1 , z> 0 
0 , z< 0 



(27) 



[0026] After judicious insertion of identity resolutions, we then obtain the results 



V[ri\ = V, [n] - ^ f J dzz" [dk [^z"e-''*^e'^y (z' )yfee'*l^*-^"l (z" \t, \k) (28) 



and 



V[n]=VAn] 



- [yzz" l7k[dz'llz"e-^'''e''^V(z')ksmik[z^^^ ^^^^ 



[0027] Next we interchange the order of the dz and dk integrals and note that 



f+oo 

dzz"i 

J-oo 



-2ikz 



" J-00 



-2/ib 



2) 



ok 



7m) 



(30) 



[0028] Consequently, Equations (28) and (29) become 

V[n\=V,[n] 

7U 



(31) 



and 



F[«]=r,[«] 



(32) 



[0029] Now the essential point to note is that when « = 0 the second terms in both Equations (31) 
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and (32) vanish since there is no derivative and k-0. Therefore, we conclude that 

F[0]=F,[0]=K.[0] 



(33) 



However, when « = 1, there is a nonzero contribution from the second term on the RHS of Equation 
(31) since one term contains "^(^) = 1 and the remaining factors are non-zero, in general. The 

OK 



second term on the RHS in Equation (32) remains zero since 



[A:sin[^(z"-z')]],.o = 0 



(34) 



and we thus conclude that 



(35) 



[0030] Thus, the approximation to the interaction produced by the Volterra formalism yields correct 
values for F[0] and V[\] while the Fredholm-based Bom series yields the correct values only for 
F[0]. hi another portion of this application, we will present a complete inversion scheme based on 
moment of Fusing the iterated Volterra-Bom series [9A]. We comment that the first order Volterra 

approximation, is not as sensitive to the near-field effects as the first order Fredholm 
approximation, Kj. 

[0031] We now turn to consider the Volterra-based inversion for interactions that can be expressed 
as sums of Dirac delta functions. 

IV. VOLTERRA-BASED INVERSE SCATTERING TREATMENT FOR SUMS OF 
DIRAC DELTA FUNCTIONS 

[0032] Rodberg and Thaler [10] have presented an interesting derivation of Fredhohn's method for 
solving the Lippmann-Schwinger equation that involves writing the interaction as a sum of Dirac 
delta functions. Essentially they argue that since 



r+oo 

F(z)= \_^dz'd{z'-z)V{z') 



(36) 



holds for reasonable interaction functions, one can write V{z) as a limit 



(37) 



[0033] While care must be exercised with such an argument (as shown in more detail in Appendix 
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B), it suggests that a useful model, especially for a scattering interaction having effective compact 
support and which is also effectively band-limited (of course, both are not rigorously possible 
simultaneously), can be taken to be 

n^)=l'j..AA^j-m^j) (38) 

where / is the finite number of delta functions needed to represent the inter-action adequately. (By 
using Hermite distributed approximating functional (HDAFs), Sj^z^z^a)) to replace the delta 
functions, we can obtain the smooth, well-behaved interaction form 



m=llSMi^-Zj\a)V(Zj) 



(39) 



about which we later make some brief comments in Appendix B. In this section of the application, 
we shall focus primarily on Equation (38).) 

[0034] We recall that the Volterra-based solution of the 1 - Z) scalar Hehnholtz equation is [7] 



^ika 

6 (•+00 



For the interaction Equation (38), this is seen to give 

e 



^ki^) = -5^+ ^Z;^.,sin[^(z^. - z)]H(zj - z)kjV{zj)i}>{Zj) (41) 

[0035] We assume, without loss of generality, that 

Zj > Zj., (42) 
In terms of the full Green's operator, we have that 

iP,{z)^{z\k)^{z\GV\k) (43) 



which, making use of the well known expansion 



0 

7=0 



(44) 



give rise to 
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v7,(z)=(z|^>+X(#on"k) 



[0036] Here, for clarity we have used ^ instead of z as the symbol to represent the ordered 
coordinate integration variables. This formula is general, but if we now introduce the interaction of 
Equation (38) the integrals can be evaluated to yield 



where through ^„ are now elements of an ordered subset of the z-points where the 5-functions of 
the interaction are located. Let N^<JhQ the number of such points greater than z, then the number 
of ordered sets of ^-points satisfying the limits on the summations is given by the binomial 

coefficiaitAy/[/i!(A^j-ii)!].Thesum X[<^, Z^^^^^ ,<f„ , which forconcisaiess we writeas Z f,,^^,...^^ 

is the sum over all N^l/[nl(N^ ~ w)!] sets. Finally, we have that 



n (z) = {z\k) + X Z (H^oH?i)(f.|^oH^2)- ■■{gn.i\GoVk){g„\k) (47a) 



11=1 



(47b) 



X sm{k[g, - z])Vig,)smik[g, - ])F(^, ). . . sin(M$-„ - g„.,])V(g„){g„\k) 



[00371 It is clear that the number of terms contributing to ip/^ (z) , for any value of z, is 



1+y^^ ^ = 2^' , which is finite (again assuming the number of delta functions in the 



interaction to be finite). Thus, for z^ Zj,N^ = 0 and only one term contributes to (z) . That is 



ipj^iz)^ {z\k) 



(48) 
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[0038] Similarly, for Zj>z^ Zj.j, iV^ = 7 and we have 

Wk (2) = {z\ k)+k^j sinik[zj - z])V(zj ){zj \ k) (49) 
and for Zj.i >z^ Zj.2, N^ = 2 and 

+ M J sm(k[zj - z])V(Zj ){zj I k) 

+ sin(A:[z^., - z])ViZj_^){zj_^\k) 

+ ^'A ,.,A , smik[Zj_, - z])V{zj_,)smik[Zj - z])ViZj){zj\k) 



etc. As one progresses in this manner from the transmission to the reflection region, the number of 
terms proUferate, but they are individually quite simple. Finally, in the reflection region the number 
of terms in the wave function is two raised to the number of scattering points in the interaction. 

Obviously, if we had such a progression of y/^ values it would be trivial to solve for the various 
V{z^ values sequentially starting from the transmission end. 

[0039] It is instructive to express this result using an A^^ + 1 dimensional vector representation where 
one component stands for the point z and the other A/^^ components represent the delta function points 
in the interaction on the trans-mission side of z. To this end we introduce a set ofN^-^\ orthogonal 
unit vectors |0}, |1}, \2}„]Nz), where |0} is the unit vector associated with z. Here we use Dirac 
notation with j/} representing the unit vector in this fmite dimensional space. Then we have that 
{s\t} = 4,/ and the identity matrix in this space is given by 

1= (51) 

[0040] We next define the matrix F by 

[j\Y\i\ = sm{k[z,-ZjW{z,)^, (52) 

Then 

(P,{z) = {0|[l+ ^112}{2|}.-[l+ ^llA^,}{iVj]|/}f \l}{z,\k) (53) 
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which provides an explicit summation of the Volterra-Bom series for this interaction. Physically, 
we see that each scattering center (i.e., delta function) in the interaction either produces a reflection 
or it has no effect. Hence the wavefunction at any point is only aware of the scattering centers that 
lie to the transmission side. The Volterra-Bom series at a point z is simply a finite sum of the 



2 possibilities. Finally, if we knew (z) in the reflection region at as many values of it as there 

are scattering points in the interaction we could in principle solve the resulting (highly non-linear) 
equations for the Vizj). We stress that such a procedure corresponds to using far-field measurements 

to determine near-field quantities exactly (essentially, one is obtaining the ) )• 

[0041] One could, of course, follow exactly the same procedure in the Fredholm case starting with 



{Pl(z)={z\k)^{z\G'v\k) 



(54) 



However, in this case there is no Heaviside-fimction in the coordinate representation of the Green's 
operator and consequently the counterpart to Equation (46) does not have an ordered sum. The result 

is that the sum of terms contributing to 1//^ (z) for any value of z is infinite. This simple interaction 
illustrates the comparative convergence properties of the Fredhohn-Bom and Volterra-Bom series. 



[0042] We next note that the T -matrix element 



(- k%\k) = (- km)= t e"''V(Zj)Ajip,(z,) 

7=1 



(55) 



From Equation (22), we have that 



F, (z) = I'jke-'^'^ t )A , iP, {zj ) 



(56) 



and thus, from Equation (50), we see that the general form of (z . ) is 

(p,{Zj)={zj\k)^R{k) 



(57) 



where the various terms in R{k) have two kinds of k dependence. First, each is proportional to a 
linear or higher power of k and, second, each contains phase factor exponentials (resuhing from 
decomposition of the sine functions). Hence, the A-integral of each term can be evaluated expHcitly. 
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The integral of the {Zj\k) terms simply reproduce the original interaction (corresponding to the first 
order Bom approximation), and, since 



/ -2ikz 



k'e 



K2 



-likz 



J 



(58) 



the terms arising fi-om R{k) contain first or higher derivatives (with respect to z) of delta fimctions. 
[0043] For example, the special case of y = 3 is 



^,(2)=lA/(z,)^(z-z,) 

+ tZ I A/,.F(z,)F(z,,)[^(2,, - z)- ^(z, 
^ y=i j'>j 



(59) 



A1A2A3 
16 



[<5"(Z3 - Z)- .y(Z2 - Z)+ ^'(Z, - Z)- J"(^, + ^3 - -^2 - ^)] 



If there are more sampling points in the interaction, the structure remains analogous but there occur 
higher order derivatives of the Dirac delta functions. Finally, we note that formally (and exactly) 

pr,+A,/2 

(60) 



Jrz,+Ay/2 
dzS'{z-z,)^0, p> 1 



[0044] It follows that averaging in the neighborhood of a sampling point averages all of the 
higher terms (/.e., the near-field or half-off-shell effects) to zero and we find that 



J»2j+A../2 
i d'zF,(z)= A,.F(z,) 



(61) 



Once the V{zj) are known, one knows the interaction. 

[0045] Of c ourse, s uch a veraging i s a f ormal e xercise f or a i nteraction t hat i s a s um o f d elta 
fimctions since both V{z) and V^(z) are not true fimctions. They only have meaning in terms of delta 
sequences. However, the operational procedure can be applied to real data to construct an HDAF- 
approximation where the HDAF is interpreted as a member of a delta sequence. Thus, the suggested 

procedure would be to take the on-shell (far-field) amplitudes (- k\T^\k) and evaluate V^(z) by 
Equation (22). Then one would use some averaging procedure such as that indicated in Equation 
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(61) to obtain approximate expressions for the A^V(Zj) on a sufficiently dense set of points to 
construct an acceptable approximation to the true interaction using Equation (39). In the process of 
this averaging, one is taking account of the effects of near-field terms in the Volterra integral 
equation. 

V. IMPLEMENTATION OF THE VOLTERRA-BASED INVERSION 

[0046] As is evident fi-om Equation (61), if we have the modified reflection coefficient (- k\T^\k) , 

we can evaluate the V{Zj) parameters approximately by a suitable averaging procedure. However, 
experiments are generally carried out under conditions that do not make direct measurement of 

k\Tf^ \kj possible. But, as shown previously, this quantity is calculated from the physical reflection 

{-klTlk) and transmission {k\T\k) amplitudes [7]. Thus to apply Equations (55), (56) and (61) 
immediately requires an additional measurement compared to a Fredhohn-based inversion. This is 
the price one pays to obtain the simphfied Volterra expressions. It is significant, nevertheless, that 
even if one cannot measure the transmission {k\I]k), the Volterra-based inversion can still be carried 
out, albeit with substantially greater required effort. This greater effort is the price one must pay to 
take account of the near-field effects in a direct fashion. 

[0047] To see how this can be done, we consider the Lippmann-Schwinger equation 



(62) 



For the interaction form of Eq.(38), this yields 



2 ^ v-y/y*v^y/ (63) 



[0048] We find the transmission coefficient 



2 , 



[0049] The reflection coefficient similarly is 
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ikk 



(65) 



[0050] The Volterra normalization is such that t. = 1 . To achieve this we note that 



h ¥k (^) = ¥l W 



[0051] Then by Equation (64), 



1 



^* ihA 



(66) 



(67) 



and 



t 



- r, 



(68) 



Thus, 



ikk 



= r. 



(69) 



and 



^ f+oo 2/ 



-2iAz 



(70) 



yielding 



V^(z)= rd(2k)e 



-2iks 2' 



(71) 



To use this expression, we substitute for (Zj ) using Equation (53), and carry out the averages 

in Equation (61) to generate a system of nonlinear algebraic equations for the (Zj ) parameters. 
The only experimental data then needed are the 
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VI. CONCLUSION 

[0052] In this portion of the application, we have considered the problem of taking account of half- 



on-shell matrix elements of either or 7]^ , We considered the moments of the interaction expressed 



in terms of and 7^ and proved. that V^[n]y /i = 0, 1 is exact, while only V^[0] is exact. This 

suggests that an inversion based on the Volterra scheme is preferred, due to the different manner in 
which the half-off-shell effects enter. This is further supported by the superior convergence 
properties of the Volterra-based inverse scattering series. We illustrated these convergence 
properties using a simple model interaction that is expressed as a sum of Dirac delta functions. It 



was shown that a formal local average of in the neighborhood of a delta function sampling point 

yields exactly the desired sampling value, AjV{zj). This is, of course, not exactly true for a real 
interaction. However, as argued by Rodberg and Thaler [IDA], using a sufficiently dense sampling, 
combined with an HDAF replacement of the Dirac delta functions, will yield a reasonable 
approximation [1 1 A]. It has the advantage that the averaging process is, in essence, taking account 
of the half-off-shell contributions since it averages them exactly to zero for any interaction 
constructed as a sum of Dirac delta functions. This delta function approach, strictly speaking, 
depends on taking a specific form for the interaction, but its application can be very general when 
the interact parameters are determined as discussed. 
A ppendix A 

[0053] The Modification of by the Factor of *^ 

[0054] Recall that the acoustic and electromagnetic wave scattering (in l-D) is of the form k^V. 
The Lippmann-Schwinger equation is 



(72) 



The transition operator is then defined as 



TM=v\i//;)=(y+vG^,k'T,)\k) 



(73) 



It is then convenient to absorb the A^-factor into the Green's s function as is done in Equation (2). 
This A^-factor has a consequence of making a Bom-expansion solution for a low-A approximation 
[6]. This is in sharp contrast to quantum scattering for which it is well known that the Bom series 
always converges for sufficiently large k [8 A, 12A]. Li the Volterra case, this A^-factor cannot 
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prevent convergence no matter how large k is. 
Appendix B 

[0055] Representation of V{z) using Dirac Delta Functions 

[0056] In a series of studies [11 A] it has been shown that a well behaved function can be 
represented to controllable accuracy by an HDAF approximation. From this point of view, the model 
interaction form 

V{z)-Y.SMi^-^jW)Vi.^j) (74) 

(see Equation (35)), with a suitable choice of the imbedded parameters K(Zy), can be made to 
approximate any realistic interaction closely. (Of course, a limitless number of delta-function 
approximations could be utilized, but as we later discuss there are advantages to the HDAF 
approximation.) If V is effectively compact (and band-Umited) then only a finite number of K(zy) 
parameters are needed. The Bom expansion of the wavefunction is given by 



ni^)={z\k)^Y.{^YG,vr\k 

n=\ 



(75) 



where y/t{z) and Go can represent either y/l (z) and or (z) and Go . Since Visa local 



operator 



{z\{G,Vy\k) = gz"(z|Go|z")K(z")(z"|(GoF)"-'|^) (75) 



and hence for this model interaction we have 



(77) 



[0057] Now SJ^z -Zj\ cr), as a member of a delta sequence, is highly localized around Zf In fact, the 
degree of localization can be controlled by our choice of the width parameter, cr. The flexibility of 
the interaction form, V{z\ of Equation (74) to represent any realistic interaction depends on its being 
"well-tempered,*' as has been discussed elsewhere [11]. This essentially requires that as a is 
decreased the number of Zy-points must be increased. By suitably localizing SfJ^z*' -zj^d) relative to 
how the rapidly the rest of the integrand varies, we can write 
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I 



[dz"{z\G, \z")s^{z"- zj I C7)(z" l(GoF)- 1^) 
= [^"(z|Go|z, [z"- Zj I C7){zj ICGoF)"-' 1^) 



(78) 



to arbitrary accuracy while still keeping the number of z^-points in Equation (77) finite. 
Since 



[dz"S^{z"-Zj\a) = 



1 



(79) 



we then have that 



(80) 



[0058] This is clearly tantamount to making a delta function approximation for V{z). Repeating this 
process for {zj\(GQV)"'^\k) etc., we can ultimately reduce the Bom expansion of the wavefunction to 
a form equivalent to that obtained starting with a interaction that is a sum of delta-functions. In this 
sense the delta-function interaction is quite general, but its application (i,e,, in the context of the 
above discussion, fixing u and the number of sample points) sensitively depends on the particular 
problem. From the basic theory of Fourier transforms, we know that the variation of the integrand 
in Equation (78) is controlled by its band-width in Fourier transform space, and, in turn, the band- 
width controls the width of the spacing between z^-points through the Nyquist relationship. Clearly 
one wants the band width to be as small as can be reasonably chosen in order to make the sample 
point spacing as large possible. 

[0059] At the same time the above discussion, also makes clear that one also wants to minimize the 
spread of the delta-function approximation about the point z, in Equation (78). The problem of 
simultaneously minimizing the band-width and the spread of the delta-function approximation in 
z-space is addressed through the choice of the delta-function approximation of the form used in 
Equation (74). It has been shown elsewhere that, in the sense of the Heisenberg Uncertainty 
Principle, the HDAF approximation (i.e., Sj^z-z^a)) is the best one can do [13A]. 
[0060] The above discussion examines the conditions under which a interaction that is a sum of 
delta-functions can adequately represent the true interaction. However, it does not justify the 
procedure of Equation (61) per se since this averaging procedure, which is strictly correct for a 
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delta-function interaction, depends explicitly on the fact that the delta function interaction is not well 
tempered. Precisely how this procedure should be implemented/modified for a reaUstic interaction 
is a subject for further research. 
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INVERSE SCATTERING THEORY: RENORMALIZATION 
OF THE LIPPMANN-SCHWINGER EQUATION FOR 
QUANTUM ELASTIC SCATTERING WITH SPHERICAL SYMMETRY 

VII. INTRODUCTION 

[0061] There have been several general approaches to the inverse scattering problem in quantum 
mechanics. The earUer of those was pioneered by Jost and Kohn[ IB] and Moses[2B] and it is based 
on the Bom-Neumann perturbation expansion of the Lippmann-Schwinger integral equation 
describing quantum scattering. Additional work on the approach includes that of Razavey[3B], 
Prosser[4B], and most importantly, in the context of the seismic inverse problem, by Weglein and 
co-workers [5B]. The key mathematical issue in the approach concerns the convergence of the 
resulting inverse scattering series, and this can be deferred, at least for some aspects of the problem. 
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by considering certain subseries[5B]. 

[0062] The other general approach has been that pursued e.g., by Marchenko[6B] and R.G. 
Newton[7B]. hi these approaches, alternative integral equations (of the Volterra-type) are derived 
leading to extremely robust behavior under iteration, Le., absolute convergence independent of 
interaction strength. So far as we can tell, the principle difficulty associated with these approaches 
is in the nature of the input data required for their implementation. Lideed, it is true in general for 
quantum scattering that experiment does not readily provide the quantities that are directly involved 
in the inversion formulae[7B]. This is in part a consequence of the fact that in quantum mechanics, 
probabiUties rather than ampUtudes are observed, thereby leading to ambiguities in phases. The 
present application is not primarily directed at dealing with this issue, although our results are 
interesting from this aspect. We shall assume that either measurements of angular distributions are 
available experimentally since these do provide the sort of phase information that one desires, or that 
one has access to integral cross sections for a range of coUision energies. 
[0063] The approach which we shall pursue has its origin in the first class of methods [1B-5B]. 
These methods are most simply formulated in terms of the solution, by iteration, of the Lippmann- 
Schwinger equation for the transition amplitude. Thus, for structureless particle scattering, in 3-D, 
one has 



T=V+VG^T (81) 



where 



Go = I > (82) 
E - T - IS 

H=T+V (83) 
denotes the non-interacting Greens fimction and the Hamiltonian, H, is the sum of J, the kinetic 
energy, and F, the interaction responsible for the scattering. We view Equation (81) now as an 
integral equation for F (rather than for T): 

V^T-VG;T (84) 

[0064] Then a power series solution for V in terms of T has the form 



00 



V=T- TG*T+ TGlTG*T-'-- = ^ Vj 



(85) 
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Such an expression is problematic since it requires knowledge of the off-shell r-matrix elements 
(which are generally not available since they are equivalent to near-field measurements of the 
wavefimction). However, for the case of a local potential, we must interpret Equation (85) as a sum 
of local, effective interactions which (provided the series converges) add up to the t rue, local 
interaction. Thus, consider the first order term: 

V, = T (86) 
an arbitrary off-shell matrix element of this then is of the form 

{k'\V,\k)={k'\T\k) (87) 

= -^l^e-^''%ine-''' (88) 



= V^(k-k') (89) 

[0065] Thus, due to the local character assumed for F(and therefore also for = 1, 2, ...), we can 
obtain all needed J-matrix elements for Equation (85) once Vy (r) is determined [1B-5B]. This 

results from the inverse Fourier transform of Equation (89), where in particular, we consider 

backward scattered amplitudes[7B]. Then -k and 

P;(2^)=(-^|71^) (90) 
V, (r ) = 2 Idke-""''-'?, (2k) (91) 
Notice that all of the matrix elements of T can be gotten firom the (2k) , simply fi-om the 



condition 



{k'\T\k") = V,C2k) (92) 



where 



k - 2 (k'-k') ^g2) 
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[0066] Thus, Equation (85) can also be expressed in the form 

V=V,- V,G*,V, + V,G*,V,G;V,— (94) 

This is all well and good except that it can only lead to well-defined results if the perturbation 
expansion Equation (89) converges. Unfortunately, this is extremely difficult to ascertain in general 
and it depends on the strength of interaction, F, the existence of bound states in the spectrum ofH 
and the energy of the collision process, etc, [7B]. In general, the expansion does not converge if the 
interaction is too strong (or if it supports bound states). 

[0067] The goal of the present application is to provide an alternative inverse scattering series 
approach which is guaranteed to converge absolutely, independent of the strength of the interaction. 
The application is organized as follows. In the next section, we derive a renormalized inverse 
scattering series and discuss its convergence. In section IX, we discuss the information required to 
apply the new inversion and in section X we illustrate the approach by applying it to a simple model 
scattering system. In section XI, we discuss our results. 

Vm. RENORMALIZATION OF THE LIPPMANN-SCHWINGER EQUATION 
[0068] We begin by remembering that Equations (8 1)-(94) also apply in an appropriately modified 
form if one considers the various partial wave components. For simpUcity, we shall restrict 
ourselves to spherically symmetric interactions in this application, but the method is general [8B]. 
The radial Lippmann-Schwinger equation is well known to be [9B] 

where !•> (tj is the usual greater (lesser) of the pair (r,r'), j, is the regular spherical Bessel 
function, is the f"* spherical Hankel function with outgoing wave condition, ^^^^ is the partial 
wave component of the scattering boundary condition solution to the Schrodinger equation, and 
G/Qj^ is defined by the second equality in Equation (95), Specifically, 

= - ^ K' tttJ w\ ) (96) 
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t 



'■ 



/ m 



(97) 



1 e 



iK\r-r 



2mk 



r-L I Yrir)Yr{ryh,{kr^)j,{krJ 



n 



I m 



(98) 



I m 



(99) 



The asymptotic form of (r) (for any r > r„„ such that V{r) = 0) is 



where 



(100) 



(101) 



With this definition, the (unitary) iS-matrix, S,. satisfies 



5/ = 1 + 2zT/'^ 



implying that 



r/'^ = e"" sin 77, 



Here rji is the usual phase shift, and 



St = e"" 



It will also prove necessary to define an additional quantity: 



2mk r ^ 



(102) 



(103) 



(104) 



(105) 



[0069] We note that in general, T/^^ is not a directly measured quantity nor is it immediately 



obtainable from measured quantities. Finally, the differential scattering ampUtude,y(6!). is 

/(^)=^Z(2/+l)^/(cos^)r, 



(106) 
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where 6 is the angle between the incident relative momentum vector, k , and the direction of 

observation, r . For mathematical simplicity, we shall also assume that the interaction has "compact 
support*', /.e., it is zero outside the range r^^x' 

^W=0, r>r_ (107) 

In general, however, our results will hold for interactions that are not too singular at r = 0 and that 
tend to zero faster than as r Following Sams and Kouri[10B] and Kouri and Vijay[l IB], 
we rewrite Equation (95) as 

2mk 

h 



1 - ^ l-A''' KMHr')y^!,(r1 



2mk 



(108) 



But 



hj (kr' ) = Tj, (kr) + ij) (kr) (109) 



so we write Equation (108) as 



2mk r , (110) 



h 



[0070] We recognize that the factor, [1+ J/^^ + /T/'^], while unknown, is simply a constant 



normalization, so that 



Kir)=u,,ir)[l+T,''' + iT,^'^] (HI) 



where 



= J,ifcr)+ Idr'r" G„,(ry)V(r')u^ir') (112) 



2f72k 

G,oM(<nr')= - —r[r?i(kr)j,ih'')- j,(h-)Tj,(kr')] (113) 
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Equation (112) for «ft(r) has the tremendous virtue, compared to the Lippmami-Schwinger equation 
for 

¥ik (^) > being a Volterra integral equation and under iteration it converges absolutely and 

uniformly for all appropriately measurable interactions. This is because the kernel, G^^j^ (r,r' ) , 

is triangular, implying that the Fredholm determinant is identically one[7B]. Consequently, it has 
no zeros and the Fredholm solution encounters no singular points. This is the best possible 
mathematical situation one can ever have! 

[0071] However, we still must address the problem of how to make use of Equation (111), since 
T^^^ is not readily available. Before dealing with this, we note that in analogy with our earlier work 



on acoustic scattering[l IB], we can introduce a partial wave transition operator, T^: 



(114) 



t, = v^vg;,j, 



(115) 



and the Volterra-based auxiliary operators: 



(116) 



or in the coordinate representation, 



2mk 



(117) 



Then we define 7] such that 



(118) 



(119) 



We see that 



r 



T ^ T 



1- 



v 



2mk 



ji 



J 



(120) 



= 7; 



1- 



Imk 



J I 



2mk 



Jl 



(121) 
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'I 



It follows that 



(122) 



We then see that 



f (I) 

^1 ~ 



l_7j(2)_ (123) 



This relation enables us to express the perturbation expansion of ^in terms of J/'^ and TP^ ,and 



then ultimately in terms of 7]^ We stress, however, that from Equation (119), 



V=T,-VG„,T, (124) 



= T, - T,G„,T, + T,G,,J,G,,J, + '- (125) 

which expression converges absolutely and uniformly independent of the strength of the interaction. 
We shall again restrict ourselves to local (and for this application, spherically symmetric) 
interactions, F(r). 

[0072] now consider how to determine the V. , defined by 



00 



7=1 



(126) 



Vi = T, (127) 
^2 = -y^G„,V, (128) 

^3 = ^i^^m^iG^mf". = -^2^?m^. (129) 



= VjA,V, (130) 



Fj'> - - ^ 1 drr'jj ikr)V, (r) (131) 
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2mk 



(132) 



We have written the upper limit as <» rather than tj^ 'm anticipation of the fact that, provided (r) 
tends to zero faster than l/i^, the more general result holds. It is not difficult then to show that 



(1) 



1/ 



1- V^^^ - iV 



(1) 



(133) 



(1) 



21 



1- V^^) _ ivO) 



(134) 



These results have some extremely interesting features. First, note that both V^f^ and V^f^ are 
purely real. Consequently, Equation (133) guarantees satisfaction of the optical theorem, since 



(135) 



(This also implies that one cannot eliminate V^P in favor of F,}^^ by using the real and imaginary 

parts of Equation (133)). Second, Equation (1 33) is sufficient to enable an inversion, provided that T/^^ 
is known. Thus, Equation (133) is viewed as an equation that is satisfied by the first order radial 



function, K (r) : 



2mk 



2mk 



(136) 



r Jo 



[0073] Now, of course, the inversion is not so simple since one no longer has Cartesian Fourier 
transforms to invert. One procedure is to express V^(r) in some basis set (e.g., Bessel functions); 

an ahemative is to use the distributed approximating functionals. In any case, it may be necessary 
to obtain (r) on a numerical grid. We conclude that the fundamental results are Equations (133)- 



(134), along with: 
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(137) 



We now turn to discuss the implementation of this approach in terms of measurable quantities. 

IX. EXPERIMENTAL DATA REQUIREMENTS FOR IMPLEMENTATION OF THE 
VOLTERRA-BASED INVERSE SERIES 

[0074] hi quantum mechanical elastic scattering, the optimum measurements are the differential 

angular distributions, which are determined by Due to azimuthal symmetry, we, in fact, 

consider 2j^0)\hin&: 



da{6) . In 

= sin^— T 

dQ k' 



X(2/+l)P;(cos^)r/^^ 



(138) 



Clearly, if sufficient number of scattering angles are measured, one can (in principle) determine the 



,An altemative that avoids having to determine the individual 7)^^^ is to use Equation (136) 
directly, along with an appropriate representation of (r) to derive a system of inhomogeneous 
(non-linear) algebraic equations which can be solved. For example, if we expand V^(r) in a basis, 



(139) 



it can then be seen that 



— — ;— = Sin 



1 ~ n'^lp'^lp' 



(140) 



where 



2mk f«> 2 , 

^ip = - -p- 1 drr j, {kr)^p{r) (141) 



and 



H.. = - 



2mk 



IP ^2 



[ drr''h;ikr)j\{kr)^^(r) (142) 
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In general. Equation (140) would be solved by a least square method, using more lvalues than the 
number of terms in the expansion over p, Equation (139). The redundancy is useful for averaging 
out noise. Another alternative is to use angular measurements at a small number of ^'s, but for a 



range of collision energies (E = h k 1 2m) to obtain an over-determined set of simultaneous 
nonlinear algebraic equations to solve. 

[00751 Another interesting approach can be based on integral cross section measurements. Thus, the 
integral cross section at energy E is well-known to be 



47r 



^(^) = TrI(2/+i)|7; 



(143) 



leading to the expression 



^(^)=-^Z(2/+l) 



^"2-1 p^xp'^ip' 



(144) 



One must evaluate a(E) at enough energies to generate the requisite algebraic equations for the V^^ . 

In expression (140) and (144), it is clear that lower energy measurements will be numerically less 
complicated because the partial wave expansion will converge with fewer angular momentum states. 
However, we also expect that (at least for potentials with a repulsive core) the short range part of 
the potential will be less accurate than the longer range, if one uses low energy data. 
[0076] It is important to note that in the case of the integral cross section approach, one makes no 
use of phase-dependent effects; indeed, only the \T^^ enter the expression. Neither the basis set nor 
D AF approaches necessarily require knowledge of the phase of the J,, but one speculates that an 
inversion based on angular measurements will be more robust (in terms of accuracy) than one based 
on integral cross section measurements. This remains to be tested. If either the integral or differential 
cross section approach proves to be feasible, this will represent an extremely attractive feature 
compared to approaches that require determination of the individual partial wave phase-shifts. 

[0077] Finally, we point out that once V^{r) is known, all higher order Vj{r) can be computed 

in the coordinate representation, using Equation (137). We notice that the same radial functions, 
Vj (r) , result no matter which partial wave is considered. This is a consequence of the assumed 
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spherical symmetry of the original potential. This provides an intemal consistency condition that 
must be satisfied. 

[0078] Of course, the above ideas, while formally correct, still must be tested on actual experimental 
data. An important issue is that of the effects of noise and inaccuracies in the data. In this regard, 
DAF-based methods may offer advantages. In a DAF approach, an approximation to the identity, 
of the form 



/2=0 



r — T 



\ a ) 



(145) 



is employed. Generally, it has the properties that 



lim tJ^(r-rV)= ^(^-^*) 



(146) 



and 



lim^^(r- r'\(j) - 5{r - r') 

<T->0 



(147) 



here o is a length-scale parameter. For finite M or nonzero C7, - r\d) is not a projector onto an 
orthogonal subspace of Hilbert space (though it is, usually, such a projector onto the Schwartz 
subspace). The DAF approximation (which can be made arbitrarily accurate) to the potential is then 



K,(r)= prV'2^^(r-r'|<7)^^(^') 



(148) 



M 



\ (J J 



(149) 



[0079] If we evaluate the integral by quadrature, then 



(150) 



and the F,^ in this approach are seen to be the discrete, point-wise samples of the furst order 

potential. Clearly, one will obtain expressions for do/dO and o{E) that are entirely analogous to 
those resulting fi-om the basis set expansion approach. Both methods require the solution of 
nonlinear algebraic equations and the same experimental data is employed. 
X. A SIMPLE EXAMPLE 
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[0080] By far the simplest scattering problem to solve from the point of view of the Lippmann- 
Schwinger equation is for a local, Dirac delta function potential. In our first study of inverse acoustic 
scattering, we found that the Volterra-based series converged to the exact result in a single term. A 
simple 3-D analogue is the spherically symmetric potential 

Vir)= JiSir-r,) (151) 
It can then be shown that the exact transition amplitude is 

-"^^0 7/(^0) 

^^'^ ^ 2mk Imic ^^^^^ 

1+ -p-'^o^^/(^o)7/(^o)+ -^^oV/(^o) 



Given such detailed input, we can use Equation (136) directly: 



2mk f« o ~ 2mk r® 

1+ dn''rj,(kr)V,(r)j,(kr) + —rldrT'jf{kr)V,(r) 



Obviously, the solution is (independent of the partial wave considered) 



(154) 



[0081] Next, we must evaluate the higher order terms in the expansion of V(r) in terms of the 
Vjir). By Equation (137), 



= -A4-ro)pr'G,o,(r,r')A4'-ro) 



(155) 



(156) 



(157) 



Clearly, due to the behavior of G/^^ (/*o .'b) > ^2 (f) identically zero, no matter what that value 
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of /. Further, by Equation (137), all higher Vj's are also zero. We conclude that 



^(0=i;^W- A^C^-A-o) (158) 



Thus, the Volterra-based inverse series again converges to the exact result in a single term. 
[0082] Of course, in general, one does not know the individual T^s. In this model problem, the 
differential scattering amplitude is 



2mkr, 



Y,(2l+V)P,(cos0)jf(kr,) 



0 / 



2mk , , (159) 
l+^^o'V(^o)y/(^o) 



and the cross section is the square of its modulus. The convergence of this partial wave series results 
from the property of the Bessel functions, j/Jcr^), that j^kr^) - 0 for / > kr^. Again, one can in 
principle obtain sufficient equations and obtain the exact result. The basic conclusion is the same, 
namely that the Volterra inversion converges to the exact result in a single term. 
XI. DISCUSSION OF RESULTS 

[0083] In this application we have presented a new approach to the inverse scattering problem in 
quantum mechanics. Although attention was focused on purely elastic scattering by spherically 
symmetric potential, the method is quite general. Indeed, it not only can be applied to quantum 
scattering, but to many other types of processes. Any process that can be described by a Lippmann- 
Schwinger type, causal (or anticausal) integral equation should be amenable to the approach. The 
method is based on a renormalization transformation of the Lippmann-Schwinger Fredhohn 
equation to obtain a Volterra integral equation. In quantum scattering, such equations are well 
known but principally used to analyze the analytic structure of the ^-matrix. An exception is the 
earlier work of Sams and Kouri[10B], who utilized the renormahzation point of view to develop a 
noniterative numerical method for directly solving for the coordinate representation of the T-matrix. 
The principal benefit of the renormalization to a Volterra equation for inverse scattering is the fact 
that their noniterative solutions converge absolutely and (under relatively mild conditions) 
uniformally independent of the strength of the interaction. This feature allows us to utilize the 
Volterra equations in a manner sunilar to that pioneered by Jost and Kohn[ IB], Moses [2B] and most 
recently by Weglein[5B], but with the guarantee that the inverse series always converges. 
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[0084] In the case of quantum scattering in 3-D, the results are complicated by the facts that (a) the 
renormalization factor is no longer a directly measurable quantity as it is for acoustic scattering in 
1-D, (b) the different partial waves do not separate in a simple fashion, (c) the equations which one 
must solve to determine the potential are nonlinear, due to the intrinsic nature of quantum 
mechanics. However, there are no difficulties in principle with the present method. Furthermore, the 
present inverse series does not require the determination of phases. It can, at least in principle, be 
applied either to differential or integral cross section measurements. If it indeed is the case that 
sufficiently accurate results can be obtained without requiring determination of phase sensitive 
quantities, this will provide a major advantage over other inversion equations for quantum 
scattering. 

[0085] For the case of scattering by a spherically symmetric Dirac delta function potential, the 
convergence to the exact resuU is obtained with a single term. By contrast, the Bom-Neumann 
inverse series based on the Lippmann-Schwinger equation yields the result 



for the first order, effective local interaction. It is immediate that any real obtained fi-om the 
above will introduce unphysical b ehavior since the left hand side of the equation is real; i.e.. 
Equation (160) manifestly violates the optical theorem for real F,(r). Comparing this to Equation 
(132) and using (152) leads to 




(160) 




2mk 



(161) 



1 + 



n 



a solution of this equation is seen to be 



F, (r) = 



(162) 



The second order correction is 



SpecincadoD 

Express Mail No.: EV 328 518 796 US 



Page 35 



ROBERT W. STROZIER. P.l_L.C. 



2mk ^ 
1+ ^"^'-^ 



^oV(^o)y/(^o) 



(163) 



^ 0 (164) 

Thus, in this case, the first order term of the series does not yield the exact answer in general, and 
it does not consist of a single nonzero term. In fact, one can then see that one must sum the infinite 
series analytically in order to obtain the correct result for values of X that are outside the 
convergence limit of the series. Equation (1 60) corresponds to the first term in the Taylor expansion 
of the denominator on the right hand side of Equation (161), which is analogous to the situation we 
encountered in our previous work on ID inverse acoustic scattering. Such an expansion converges 
only for sufficiently small -values (as well as also depending on the value of r^. Of course, it does 
permit one to sum the infinite series analytically to obtain the result that holds outside the 
convergence limits of the series itself[l IB]. As is also usual for the Bom-Neumann expansion in 
quantum scattering, the approximation does eventually converge for high enough energy, E (large 
enough k). 

[0086] We are currently exploring the inversion of quantum 3D elastic scattering by a nonspherical 
target, as well as various other wave phenomena. Of particular interest are the cases of acoustic and 
electromagnetic scattering in fiiU 3D. In addition, we shall carry out test calculations to verify that 
one can use non-phase sensitive, integral cross sections to carry out an inversion. These results will 
be reported as they are obtained. 
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INVERSE SCATTERING THEORY: RENORMALIZATION OF THE 

LIPPMANN-SCHWINGER EQUATION FOR 
ACOUSTIC SCATTERING IN ONE DIMENSION 

XII. INTRODUCTION 

[0087] The inverse scattering problem has enormous importance both for practical and theoretical 
applications. The former include hydrocarbon exploration and production, medial imaging of many 
varieties, nondestructive testing, target identification and location, etc. The latter include relating 
interactions governing atomic and molecular systems to experimental measurements, determination 
of the structure of surfaces and condensed matter systems, imaging of nanostructures, etc. In much 
of the literature, the focus has been on determining the conditions under which the data inversion 
will yield a unique result and precisely what information is required to make an inversion possible. 
In terms of algorithms employed for various types of imaging, an important practical tool is the first 
Bom approximation, which assumes that all scattering is direct, involving a single interaction of the 
probe with the target. Of course, this is known to be incorrect. Indeed, most imaging procedures or 
algorithms typically make use of some assumed model for the propagation of the probe signal or 
disturbance in the scattering medium. 

[0088] Generally, inversion is practical only in their circumstance that there is a sufficiently small 
difference between the propagation of the probe signal within the target and its "reference 
propagation" (low contrast between the target and the reference medium). Over the last decade, 
Weglein and co-workers [ 1 C] have pioneered inverse acoustic scattering methods that do not require 
an assumed propagation velocity model within the medium. Their approach is based on the early 
work of Jost and Kohn [20], Moses [30] and Razavy [4C] that used the Bom-Neumann power series 
solution of the acoustic Lippmann-Schwinger equation, and a concomitant expansion of the 
interaction in "orders-of-the-data". Reversion of the Bom-Neumaim series leads to an order-by-order 
scheme for evaluating the terms of the series representation of the scattering interaction in terms of 
the measured data; e.g., only the on-shell reflection amplitude is required to invert for a local 
interaction. In principle, the method is completely general and requires no prior information about 
the target or the propagation details of the probe signal within the target. 
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[0089] The only fundamental limitation of the approach appears to be the finite radius of 
convergence of the Bom-Neumann series solution of the acoustic Lippmann-Schwinger equation. 
This is generally analyzed using the "spectral radius" of the Fredholm kernel of this equation [Morse 
and Feshbah, [5C]; Newton, [6C]], and in particular by theL^-norm of this kernel. References and 
very clear discussions of the issues involved in the convergence of the Bom-Neumann forward 
scattering series can be found in [Goldberger and Watson [7C]; Newton [6C]]. Despite this 
limitation, Weglein and co-workers [1] have made significant progress using this approach by 
introducing the idea of "subseries" within the Bora-Neumann expansion, which are associated with 
specific inversion tasks. This expresses the inversion series in terms of a set of subtasks which can 
be carried out separately firom one another. 

[0090] A particularly significant benefit of this approach is the fat that the convergence properties 
of the subseries studied to date are much more favorable than those of the fiiU Bom-Neumann series. 
Indeed, empirical evidence has been very encouraging regarding the convergence of the inverse 
series. However, the nature of the kernel of the L ippmann-Schwinger e quation, viewed as an 
equation for the interaction in terms of the J-operator, is such that its maximum eigenvalue always 
depends on the explicit nature of the on- and off-shell J-matrix and general statements regarding 
convergence are difficult to obtain [Prosser [8C]]. 

[0091] Another, more robust approach to solving integral equations is that due to Fredhohn [9C], 
which can be viewed as a generalization of the well-known Cramer's method for solving systems 
of linear simultaneous algebraic equations. Consequently, fimdamental to the approach is a 
continuous generalization of the determinant of coefficients and its minors. Under the circumstances 
that the integral equation is of the Volterra type, the "Fredholm determinant" can be shown to equal 
one and the Fredholm solution reduces to a Bom-Neumann expansion, all-be-it one that converges 
absolutely independent of the scattering interaction strength. Consequentiy, for such VoUerra 
equations, the Bom-Neumann expansion possesses the most robust convergence properties for 
which one can hope. 

[0092] Some years ago, Sams and Kouri [IOC] (for noniterative computations in quantum 
scattering) and Kouri [IIC] (for electromagnetic scattering) showed that one could carry out a 
renormalization transformation of the Lippmann-Schwinger equation into a Volterra equation form. 
Although the Volterra equations for quantum scattering were well known [Goldberger and Watson 
[7C]; Newton [6C]] , previous studies had focused almost exclusively on their use for studying the 
analytic stmcture of the iS-matrix and the scattering state. The work of Kouri and co-workers 
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concentrated on making use of the Volterra form of the scattering equations to create a noniterative 
computational algorithm. Their approach, however, made essential use of the "triangular" character 
of the Volterra equation kemel, which in one dimension (ID) is 

^(z z') = 0, z ^ or K(z, z') = 0; z ^ z' (165) 
combined with a Newton-Cotes quadrature to solve the equations by a noniterative recursion. 
However, it is also well-known that the property, Equation (165), underlies the extremely robust 
nature of the convergence of these Volterra equations with respect to an interative solution [Morse 
and Feshbah [5C] ; Newton [6C]] . Indeed, the Bom-Neumann series solution of the Volterra equation 
converges absolutely, irrespective of the magnitude of the (in general complex) coupling strength 
of the interaction! Furthermore, the convergence depends on the global behavior of the interaction 
(essentially whether it is measurable in a particular sense) and not on its smoothness. For ID 
interactions having compact support (and for even more general interactions in the case of 3D 
scattering), the iterative solution of the Volterra equation converges uniformly on any closed domain 
of definition in the scattering position variable. Again, under certain relatively weak conditions on 
the interaction, the iterative solution is an entire function of the scattering wave nimiber, k [Newton 
[6C]]. 

[0093] Thus, the possible benefits of formulating acoustic scattering in terms of Volterra kernels 
appear substantial. The infinitely large radius of convergence of the Bom-Neumann series solution 
of the Volterra equation is of especial interest fi-om the standpoint of the inverse acoustic scattering 
approach o f Weglein and co-workers [ IC]. It seems natural, therefore, to i nvestigate p ossible 
benefits of using the renormaUzation technique as a firamework for developing an inverse scattering 
series. In fact, we shall show that it is possible to establish general, rigorous convergence properties 
for the inverse acoustic scattering series for the first time, and in the process show that its radius of 
convergence is also infinite! We shall restrict our discussion here to ID scattering but our approach 

ft 

is completely general and extends to higher dimensions [12C]. 

[0094] This portion of this application is organized as follows. In Section Xm, we discuss 
renormaUzation of the Lippmann-Schwinger equation for acoustic scattering and introduce an 

auxiliary transition operator, T . This is used as the firamework to analyze the convergence of the 
forward scattering Bom-Neumaim series. The approach is illustrated by applying it to scattering by 
a Dirac delta fiinction model interaction. In Section XIV, we show the relationship between the 

interaction as a fimction of the physical J-operator and as a function of the auxiliary T -operator. 
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We next analyze the non-local nature of T in the coordinate representation, and then use the results 

to establish the convergence properties of the Volterra-based Bom-Neumann inverse series for the 
interaction. We include in this Section an application to the Dirac delta function interaction. Next, 
in Section XV, the Volterra inverse series is applied to the case of sound scattering by either a 
square well or barrier. Our conclusions are given in Section XVI. 
XIIL RENORMALIZATION OF THE LIPPMANN-SCHWINGER EQUATION 

A. Derivation of the Renormalization Transformation and Auxiliary Transition 
Operator T 

[0095] We assume that the reader is familiar with the acoustic scattering Lippmann-Schwinger 
equation for the transition operator, J, given b y [ Razavy [4C]; Goldberger and Watson [5C]; 
Newton [6C]] 



(166) 



where Gq^ is the causal free Green's operator, multiphed by a factor of ifc^, 



E- Hq + is 



(167) 



I^^E (i.e., k is the frequency associated with the incident acoustic wave), Hq govems the "free 
propagation" of the acoustic wave, and yV is the interaction responsible for the scattering, with y 
being the coupling parameter characterizing the strength of the interaction, hi general, /is complex. 
The additional factor of results from the fact that in acoustic scattering (as in general for 
scattering governed by a Helmholtz type wave equation), the interaction responsible for scattering 
depends on ^. The frill acoustic wave propagation (scattering process) is thus governed by the 
operator H, 



(168) 



The present ID acoustic scattering problem in the coordinate representation leads to 

r(z,z')= ;kF(z,z')+ J_ dz''yV{z,z'')\ dz'''G^,(z\z'")Tiz'",z') (169) 

By incorporating this factor of into the Green's function, we are able to treat the remaining portion 
of the interaction that depends purely on the spatial variation of the scattering interaction. Initially, 
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we restrict ourselves to "local scattering media", so that V(z,z') = V{z)Siz - z') and therefore 

r(z,z') = rf^(z)c5(z- z')+ rcfe"Go\(z,z")r(z",z') 

w—cti 



(170) 



[0096] The non-local character of the causal free Green's function, Gq^ (z,z' * ) , reflected in its not 

commuting with is responsible for the fact that T(z, z") is also generally non-local; i.e., it is 
never diagonal in the coordinate representation (except for a local, Dirac delta function interaction, 
y(z^ z') = V{z)Siz - z') = >l3iz - z')Siz - Zq)). For ID causal scattering boundary conditions, 

Gqj^ {ZyZ^ * ) is exphcitly 



Go;(z,z^')=-|e'^'--' 



(171) 



The general scattering amplitude is determined by the matrix elements of the J-operator, usually 
computed in the momentum representation, T{k\ *"), given by 

T{k\k'') = {k\T\k'') (172) 

where in general, k\ k" and the on-energy-shell wave number, k = 4e need not be equal to one 

another. The physical "reflection scattering amplitude", denoted i*(A), results when \k\ = |/r"| = |k| 
and K = -k\ 

r{k) = {rikn){-k\T\k) (173) 
In ID scattering, one can also identify the transmission amplitude, t{k\ given by 

t{k)^\^{-ikn){k\T\k) (174) 
[0097] In the work of Sams and Kouri [IOC], the renormalization transformation to a Volterra 
equation results from eliminating the \z - z"|-argument in the free Green's function in Equation 

(170). This is done by dividing the integration over z" into segments from -<» to z and from z to <»: 



T{z,z') = YV{z)5{z- z')- —yV(,z)^ dz''e"^'-'"^T(z" ,z') 



- y r^(z)f °°rfz''e-'*^^-^">r(z'\z') 



(175) 



One then adds and subtracts - ^^}Viz) ^fe''^'*^''^"^r(z'^z*), and after simple manipulation, 



one 
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« 



I 



obtains 



r(z,z') = rviz) 



5{z-z')--e 

2 



ikz 



rdz"e-'^"T(z'\z') 

J_oo 



ik 
T 



rV{z) r dz"\e-''^'-'"^ - e'*<^-^"^lr(z",z') 



(176) 



It can be verified that this is equivalent to writing Gq^ (z,z' ' ) as 



(177) 



so that 



Go,(z,z") = - - = A; sin[^ (z"-z)], z < z" 



(178) 



= 0, z^z" 



In abstract operator notation, this is 



(179) 



(180) 



[0098] This relation is extremely usefUl in oxir subsequent analysis and we shall make much use of 



it. Notice that the Green's operator Gq^ differs firom the usual causal one, Gq^ , by a solution of the 



homogeneous equation [Newton, [6C]]: 



(E-H,)G;, = k 



(181) 



(E-H,)G,, = k' 
{E- H^)[-ik7i\k){k\]=[-ik7i\k){k\{E- H.)=0 



(182) 



(183) 



The abstract version of Equation (1 76) results from substituting Equation (1 80) into Equation (1 66): 

T=yV[\-ik7^k){k\T\^yVG^J . (184) 
[00991 Next we note that the action of T on the initial state \k) is of the form 



71^) = 7V[\- ik7r{k\j\k)]k)^ jVG^.nk) 



(185) 
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(unknown) constant, c^, as 



c, = l- ik7r{k\j]k) = t{k) 



(186) 



we see that 



'l\k)=)Vc,\k)^)VG,,T\k) 



(187) 



The relationship between 71*) and the Lippmann-Schwinger pressure state, j , is 



(188) 



and thus 



(189) 



Clearly, the factor is simply a normalization constant and one can define an auxiliary pressure 
state vector \p/), in relation to ^ , according to 



\p,)=^f2^\k)^G,,)V\p,) 
[0100] The coordinate representation, {z\p,) = p/,(z) satisfies 



(190) 



(191) 



(192) 



which is recognized as an inhomogeneous Volterra integral equation of the second kind. We remark 
here that Volterra equations involving improper limits (/. , ±«>) still converge absolutely for | /|<«>, 
but they must satisfy additional restrictions on the z-dependence of the interaction. This is especially 
true in order for their iterative solutions to converge uniformly on any closed interval [zu Zj]- It is 
sufficient that the interaction V(z) have compact support and \V{z)\ be measureable. It remains true 
even for infinite ranged interactions so long as they decay sufficiently rapidly and are not too 
singular. This is discussed for similar Volterra equations in [Goldberger and Watson [5C]; Newton 
[6C]]. Throughout our discussion, we assume that such conditions are met. By Equation (190), 

results from renormalizing P^ ^ according to 
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I 



Pk 



)- 



(193) 



[0101] In fact, Ck is essentially the inverse of the Jost function [Newton [6C]]. We remark that the 
above expression also provides the physial interpretation of the "Volterra pressure wave", /;^(z) 
[13C]. Clearly, it represents a wave produed by an incident plane wave having an amplitude equal 
to 1/c^ = \lt{k). This leads to a reeted wave with the amplitude r{k)lt{k) and a transmitted wave with 
amplitude exatly equal to one. Of course, such an incident wave cannot, in general, be created 
experimentally since it requires advance knowledge of the effect of the scatterer in the fomi of 
\lt{k). However, this does not alter the interpretation of the wave Pkiz). 

[0102] Let us now retum to Equation (184), and define an auxiUiary transition operator, T, 



according to 



(194) 



It can be verified that 



(195) 



and this is the fundamental equation which will be used to analyze the inverse series for yV. (Note 
that the operator inverse, [1 - ikn\k){k\T\\ should always exist. This essentially requires that the 
operator ikn \k){k\T not have any eigenvalues equal to +1 . A worst-case would correspond to the 
inverse of T being equal to ik7i\k){k\, which cannot occur since J does not commute with J^o while 
\k){k\ does.) It is instructive to evaluate explicitly the normalization constant, c*, in terms of the 
solution of the Volterra equation. This can be quite done by combining Equations (1 85), (1 86) and 
(194) to write 



expressed as 



'I\k)= cj\k) 



(196) 



= 1- ik7r{k\T\k)c 



(197) 



so that 



1 



Cu = 



1+ ik7r{k\T\k) 



(198) 
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I 



Thus, the renormalized or auxiliary pressure state, \pi), is given by 

\p,)=M[l-Mk\T\k)] 

[0103] The physical reflection amplitude, is given by 



r(k) = -ik7r{- k\T\k)=^ -t{k)ik7r{- k\T\k) 



(199) 



(200) 



These relations provide us with the necessary tools to express auxilliary amplitudes in terms of the 
physical amplitudes. 



B. Convergence of the Born-Neumann Series for and T 
[0104] On one hand, the convergence of the Bom-Neumann series for either ) or J is well- 
known to depend critically on the size of the coupling constant, /(or equivalently, on the size of the 
"contrast" between the propagation under /Tq and that under fiT = HQ+l^yV) [Goldberger and Watson 
[5]; Newton [6]]. On the other hand, it is also well-known that iterative solutions of either Equation 

(192) or (195) converge absolutely for | /\< «> [Newton [6]]. Furthermore, the iteration of Equation 

(192) converges uniformly on any closed domain of z (for a wide class of interactions). It is useful 
to stress the origin of this robustness since it turns out to be the basis of the convergence of the 
Volterra-based inverse series for yV. The kernel of Equation (192) can be written (for all z, z") as 



= 0, z> 



(201) 
(202) 



[0105] According to the discussion in [Newton [6]; see also those in Rodberg and Thaler [14] and 
Mathews and Walker [15]] one characterizes the convergence in terms of Fredhobn's method of 
solution. This method is the continuum analogue of solving a linear system of algebraic equations, 
and it expresses the inverse of the integral kernel in terms of the ratio of the first Fredhohn minor 
to the Fredholm determinant. A, The determinant A can be expressed as an infinite series of the form 
(for the acoustic case) 



00 



(203) 



n=0 



where 



(204) 
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It is not difficult to verify that for the Volterra kernel, K{z^ z"), above. 



(205) 



and consequently, for such kernels, 



A = 1 



(206) 



regardless of the strength of the scatterer, Furthermore, by use of Hadamard's theorem [6C, 14C], 
it then can be proved that the infinite series for the first Fredholm minor converges absolutely and 
uniformly for y in the entire complex plane (more details are given in the Appendix C). It also has 
been established [Mathews and Walker [15C]] that when the Fredhohn determinant equals one, the 
Fredhokn solution is identical to the Bom-Neumann iterative solution of the integral equation. We 
conclude that iterative solutions of Volterra integral equations possess the most robust convergence 
possible. While it is true that these convergence properties are independent of the strength of/, there 
are conditions on the analytical structure allowed for the scattering interaction. These have to do 
with the integrability of any singularities and the behavior at infinity. They are discussed in [Newton 
[6C]] in some detail. If the interaction has compact support, and is not too singular, then the 
convergence is of the strongest character (z.e., absolute and imiform, leading to entire functions of 
wave number k and coupling f^. 

[0106] It is therefore clear that the essential property of the Volterra kernel is that it satisfies 
Equations (201)-(202); as noted in [Newton [6C]], this is the continuous version of the "triangular" 
property of matries. It is equivalent to the property that the Fredholm determinant of Equation (192) 
or (195) is identially one. Furthermore, it ensures that the Bom-Neumann series for /^^(z), obtained 
fi-om Equation (195), is uniformly convergent on any closed domain of z, for a wide class of 
interactions. 

[0107] We stress that this is all well-known. We have included it explicitly here because its 
implications for the inverse scattering series determining yFhave never been explicated. In addition, 
it is perhaps not appreciated that the original Lippmann-Schwinger equation itself can be directly 
iterated in a fashion that is also everywhere absolutely convergent! Obviously, such an iteration 
must differ firom the straightforward iteration of the Lippmann-Schwinger equation. 



P/ ) = yf27t\ k) + G^^)V\ ) , which leads to 




(207) 
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the proof of whose convergence depends on the Zr^-norm of the kernel, Gq^;^F . In fact, wean 
simply iterate Equation (1 89) for the Lippmann-Schwinger state ) : 



00 



00 



«=0 n=0 



(208) 



We stress that even though is unknown in Equation (208), it is simply a number and can be 
calculated directly from the known iterate- vectors, [Gli^yVX \ k) • Thus, 



1 



c = 



00 



(209) 



«=0 



This is equivalent to the iterative solution for but the point we wish to stress is that the standard, 
physical Lippmann-Schwinger equationan be iterated in an absolutely convergent fashion, 
independent of the strength of the interaction. Ofourse, this is simply a reflection of the fact that the 
Lippmann-Schwinger equation is neither purely a Fredhohn or Volterra equation. Therefore, it can 
manifest the convergence characteristics of either, depending on the manner in which it is written 
and iterated. 

C. Illustrative Example 
[0108] It is helpful toonsider an example problem in order to appreiate better the vast difference in 
convergence between the Bom-Neumann series based on the Lippmann-Schwinger Fredhohn 
equation and the renormalized Lippmann-Schwinger Volterra equation. A convenient and simple 
model scattering interaction is the Dirac delta function: 



The solution to the Lippmann-Schwinger equation is found from noting that 



(210) 



(211) 



This implies that 
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so the exact solution is 




The Bom-Neumann series solution is given by 



P; (z) = e"^ - 



iky 



/ 



1- 



iky ( iky^ ^ 



+ 



^ 2 ; 



+ 



(214) 



It is clear that the convergence of this series is determined by the requirement \k)42\ < 1, which is 
just the condition for the convergence of a power series expansion of (1 + ikyfiyK It is also evident 
that the Bom-Neumann series is convergent only at low energies in this case since k = 4e , and 
therefore only for sufficiently low E will th e convergence condition will be satisfied. This is the 
opposite of the usual situation that applies to quantum scattering [Goldberger and Watson [5C]]. Of 
course, in this simple example, one can recognize that the series can be analytically summed to yield 
the exact result valid at all k and y- In general, that will not be the case. 
[0109] We next consider the Bora-Neumann series for Equation (192) then becomes 



+ A:'J ife"sin[A:(z"-z)];'J(z"-Zo)J.. f/z'"sin[A:(z'"-z")Mz'"-Zo)e 



lit 



(215) 



We see that all terms higher than first order in the interaction vanish identically due to the 
appearance of the factor sin[Ac(zo - z^]. Thus, for z ^ z^, we obtain exactly 



ikz 



(216) 
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I 



■ 



and for z < Zo, we obtain exactly 



Pk (^) = ^'^^ + sm[k{zQ - z)\e 



ikZf 



(217) 



[0110] This does not complete the analysis since we must also evaluate using only infonnation 
generated by the iterative solution for p^{z). This is simple using Equations (196)-(200) and yields 



1 



00 



(218) 



1 



(219) 



But by Equation (216), this gives 



1 



* iky 



(220) 



Therefore the Bom-Neumann series for (z) based on the renormalized Lippmann-Schwinger 



equation results in 



,/feo 



r iky 2]' 

2 ; 



(221) 
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[0111] We conclude that the Volterra-based iteration converges to the exact answer with just the 
first order term, ail higher terms being zero. The fimdamental difference between the Volterra and 
Fredhohn iterated expressions is that the former does not involve a power series expansion of the 
normalization factor, c w hose c onvergence w ould h ave r equired t hat | kyf2\ b e 1 ess t han o ne. 

Instead, Cf^ has been factored out by renormalizing from (z) to We emphasize that this 

renormalization follows for any scattering problem that is expressible in terms of Greenes functions 
Gq^ , since it is true in general (for ID scattering) that 

Gok = G^k+ik7i\k)(k\ (223) 

We also note that analogous relationships have been derived for 3D scattering Green's functions. We 

now tum to consider the inverse scattering series for the interaction, yV. 

XIV. THE INVERSE SCATTERING SERIES FOR 

A. Fredholm and Volterra Born-Neumann Series for jV 

[0112] We begin by estabUshing that distinct Bom-Neumann series for /Fcan be obtained from 
Equation (166) for Tor Equation (195) for f . We solve Equation (166) for yFas 

rV=T{uG;,Ty' (224) 

= r(l + Go, r - ikT^ k){k\Ty' (225) 

But Equation (194), this yields 

yV = r([l - ik7i\ k){k \t] + Go, f [l - ik4 k){k |r])'' (226) 

= Til- ikMk){k\T]'\u G.jY' (227) 

so that finally, 

^= f(UGo,f)"' (228) 

It follows that, provided they converge, yV can be obtained from either of the following Bom- 
Neumann series expansions: 

00 

rf" - 1 t{- G^jY (229) 
«=o 
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00 



n=0 



(230) 



[0113] The convergence properties of Equation (229) depend, of course, on the spectral radius of 
the kernel T , which in turn depends crucially on both the on- and off-shell elements of the T- 

matrix. For this reason, general conclusions regarding the convergence of Equation (229) have been 
extremely difficult to obtain despite heroic efforts [Prosser [8C]]. We shall see that this is not the 
case for Equation (230)! 

[0114] The convergence properties of Equation (230) will be studied using the Fredhohn method 
of solvmg Equation (228). To do so requires knowledge of the properties of the kernel T , which 

are yet to be established. It is clear, however, that when both expansions converge, they must agree, 
since convergent power series yield a unique resuU [Kaplan [16C]]. In order to investigate the 



convergence of Equation (230), we now consider the non-local character of T in the coordinate 
representation. 



B. Non-local Character of T (z^z^ 
[0115] The aim of this subsection is to establish that T (z,z') = 0 when z > zMt is not difficuh to 



show that Equation (195) has the solution 



(231) 



where 



G=Gq + G^)¥G 



(232) 



(see R.G. Newton, 1982, pp. 343-344; especially Equation 12.42 and the following unnumbered 
equation). From Equation 12.40a in Newton, we see that 



(233) 



where \//*(k,z)'\s the regular (physical or causal) scattering solution of the interacting Shrodinger 
equation andy(A:,z) is an irregular solution of the same equation, introduced by Jost [6C]. Then 
defining an interacting Green's function, G(A:;z,z') , that vanishes for z ^ z', it is easy to see that 
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G{k;z,z')= G*{k;z,z')+ ky/\k,z')f{k,z) 



(234) 



[0116] Obviously, kif/^ {k,z')f{k,z) satisfies the homogeneous interacting-Green*s fiinction 



Schrodinger equation. Then 



G(z,zO= 0, z>z' 



(235) 



and for local potentials 



T{z,z')^ yV{z)5{z-z') 



+ 



t?^" \Z^z'''^{z)d{z- z'')G{z'\z''')^ ^^^^^ 



or 



r(z,z') = rV{z)d{z- z')^ y^V{z)G{z,z')V{z') 



(237) 



It is therefore clear that as a function of either zorz\ T{z^z^ ) has support determined by F(z) or 

Also, for z>z\ the first term on the R.S. of Equation (237) is zero due to the Dirac delta 
fimction and the second term is zero due to G{z,z^) . Therefore, we have proved that 



r(z,z')=0, z>z' 



(238) 



Finally, T{z,z^) has an inerrable singularity at z = z 



C. Convergence of the Inverse Series for jV 

[0117] We note next that the kemel of the Volterra-based Bom-Neumann series for )V, Equation 
(238), is given by 



Go J 



(239) 



It is necessary to compute Tr^K" ) , but it is sufficient to examine tAk^ ) to see how the general 



case behaves: 



Tr(K')= \ dz\ dz^K{z,z')K{z\z) 

•r— 00 CO 



(240) 
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This can be written as 



^K^' ) = • • • rX G,, [z, , )f[z^ , Z3 )Go, (z3 , )r(z4 , ) (24 1) 



However, by the Volterra property of Gq^ and T , a nonzero contribution can only occur if 

Zj > Z4 > Z3 > Z2 > Zj (242) 

which clearly is never satisfied. We conclude that the Volterra property is satisfied for the product of two 
(or more) Volterra kernels and 

rr(F)= 7>(^")h 0, n> 2 (243) 
It is similarly easy to prove that Tr{K) = 0 . 

[0118] We therefore conclude that the Fredholm determinant for Equation (228) equals one. This 
guarantees that, for a n ot-too-singular, 1 ocal i nteraction h aving c ompact s upport, t he V olterra-based 
inverse scattering series converges absolutely and uniformly independent of the strength of the interaction! 
This is an amazing result since it ensures that this inverse scattering series always converges for any 
magnitude (complex) coupling constant. 

D. Utilization of the Volterra Inverse Series for yVin Orders of f and the Relation to 
Data Requirements 

[0119] In order to use the new Volterra mverse series to determine yV, the final step is to develop explicit 
expressions for it in terms of "far-field" measured quantities. The standard Bom-Neumann inversion of 
the Lippmann-Schwinger-based approach, to obtain a local potential, requires knowledge only of the 
reflection amplitude, r{k\ as a fimction of k. We shall see that additional data are required in order to use 
the Volterra inverse series. Recall that by Equation (230), 



00 



t[-G^J) (244) 



/i=0 



where 



T= T[\-ik7^k){k\T] (245) 



We shall express jVasa, power series in orders of T : 
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I 



00 



n=0 



(246) 



where 



Next, recall that by Equation (195), 



00 



(247) 



(248) 



CO f CO 



n=0\ J=l 



CO 



V V, 

J y=i 



(249) 



since T is first order in /I . We then collect coefficients of each power of : 



T=V, 



0 = V + V G V 



(250) 



(251) 



(252) 



etc. Matrix elements of these expressions are first evaluated in the ^-representation and the results 
subsequently transformed to the z-representation. This is because the starting expression involves the k- 

representation matrix elements of T . However, using the lower order operators, , to express Vj solely 



in terms of and Gq^ , one can find that, in general, 



V = -V G V 



(253) 



[01 20] This is the most convenient form with which to evaluate the higher order corrections. The Volterra- 
based expressions can be compared to the Bom-Neumann inverse series based on the usual Lippmann- 
Schwinger equation [1C-4C]: 
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00 



00 



^ = I r(- Go* = Z ^% 

«=0 j=l 



(254) 



and this leads to 



00 



/7=0 



(255) 



rr. ( rr. 
00 00 



«=0V 7=1 / 



/I 



00 



/=1 



(256) 



implying then 



A': 



(257) 



(258) 



etc. Again, one can show that in general. 



V = - V G* V 



(259) 



(260) 



However, it is crucial to recognize that 



V. t V. 

J J 



(261) 



because they correspond to orders of completely different parameters ( V, is yth order in T while V is 
yth order in 7). By Equation (245) above, it is clear that each separate factor of J involves all orders of T 



and vice versa: 



[l-ik7rT\k){k\]T= T 



(262) 



so that 



T= [l-ik7rT\k){k\]-'T 



(263) 
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* 



I 



00 



= X (- ik7tf\ k){k D" f 



«=1 



(264) 



Thus it is clear that Vj and cannot be the same. 

[01211 Now we ask how can one combine measured data with the Volterra inverse series? We compute 
the back-scattering matrix element of Equation (250): 



{- k\f\k)^ t{- Kk)={- k\v,\k)^vX- k,k) 



(265) 



But T{-k,k) is not directly measured. The far-field quantities typically measured are rik) = 

-ikp(-kjt), the reflection amplitude, and tik) = \-ik7iT{k,k), the transmission amplitude. By Equation 
(245), we write 



r(- k,k) = f(- k,k)- ik7iT{- k,k)T{k,k) 



(266) 



so 



(267) 



This expression is inverse-Fourier transformed to the space-domain, yielding F, (z) . The result is 



{k) 



(268) 



One obtains the higher order according to 

F,(z)= J d{2k)e- 



- k 



V G V 



(269) 



1 f 



f>(2*)e-^ [> • e-(='«">f^., (z-) G„ (z' ,z' •)^^ (z' ') 



(270) 
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Again, it is instructive to carry out the application to scattering by the Dirac delta function interaction 
discussed earlier. In that case, we have 



(271) 



t(k) = 



(2 + iky) 



(272) 



It can then be shown that 



(273) 



[0122] The second order corrections is given by 



Kiz) = f>(2^)e-^''^ Idz^rdz-e^^'^'^^^^^^^ (274) 



= r 



[742^)e"^'^Go,(zo,Zo)-0 



(275) 



It should be clear that all Vj vanish for j ^ 2. We conclude that the Volterra inverse scattering series 

converges to the exact answer in a single term, in the same manner as the forward Volterra series for the 
Dirac delta function interaction. We note that a crucial change from the Bom-Neumann approach to the 
Lippmann-Schwinger-based inversion is that now, we require both r{k) and ^(A:) to use the Volterra 
inverse series! 

[0123] Before leaving this example, we point out that Razavy [4C] has considered the interaction 2XS(x) 
within the Lippmann-Schwinger-based inverse scattering series. Again, analytical results are obtained 
and the exact result is equal to the sum of the Fj and Vj terms. However, the higher order terms were not 
evaluated, they do not vanish, and involve alternating signs. Thus, it appears that the series is only 
conditionally convergent, depending on the order in which the terms are grouped and summed. 

E. The Volterra Series for Non-Local Potentials 
[0124] Up to this point, we have established that the Volterra property is shared by the coordinate 

representation matrix element, (^z\t z') , if the interaction is local. In fact, we now show that this result 
is true for all non-local interactions so long as they also possess the Volterra kemel property. We recall 
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that the exact solution for T is 

f^yV^-y^VGV (276) 

where 

l-W 

G = Go + Go;^G (277) 

and 

G{z,z')= 0, z> z' (278) 

Now we do not restrict Fto be local but we require that 

F(z,z')=0, Z>z' (279) 

[0125] We want to prove that it remains true that r(z,z') = 0, z > z' . Our equation now is 

f(z,z')= r^(z,z')+;^' r^-f dz''T{z,z'')Giz'\z''')V{z''\z') (280) 

where we use the facts that ^" must be greater than z and must be greater than z'", due to the presence 
of the factors V{z, z") and V{z''\ z"). Now suppose that z>z'. The first term on the RHS of Equation 
(280) vanishes for this condition. But the second term only has non zero contributions for > z''^ > z'' 

> z since all terms involving z " > z vanish due to the factor G(z' ' , * ') . Therefore, the only nonzero 

terms contradict the condition z > z'. We conclude that /'(z, z' ) = 0 if z >z\ if the non-local potential 
has the same property. 

[0126] The analysis that the kernel of the inverse scattering equation has the Volterra property and, 
therefore, a Fredholm determinant that equals one can be carried out and we do not write it explicitly 
here. If one has the most general form of non-local potential, F(z, z\ for ID scattering, then it tums out 
that the inversion requires measuring both the far and near fields. Clearly, the operator equations obtained 

from the inversion of ^/Kin terms of the Vj 's hold regardless of whether the potential is local or non- 
local. This is also true of the absolute convergence of the series for k^yV in terms of the Vj (provided V 
itself has a Volterra-kemel structure). However, if the non-local potential has the general form F(z,zO, 
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then Equation (250) becomes 



f{k\k)=vXk\k) (281) 



with A' and A independent of one another. Then one uses Equation (245) to determine the off-shell 



elements of T in temis of the physical J-matrix elements: 



T(k' , k) = f{k' , k) - iknf{k' , k)T{k, k) (282) 



so that 



\-ik7tT{k,k) (283) 

Thus, knowledge of the on- and half-off-shell 7-matrix elements enables one to determine the 
corresponding elements of the T -matrix. Then inverse Fourier transforming on both k' and k 

independently yields Kj (z,z') , which enables one to determine all higher Vjiz^z') J>\. 

[0127] It should be clear that scattering interactions that can be expressed in the form 

J" 

V{z) 5{z - z' ) will also produce Volterra kernels (that is, interactions involving derivatives of the 

field). Thus, the range of systems for which our results hold is very broad. 
XV. APPLICATION TO THE SQUARE WELL OR BARRIER 

[0128] As a second example, we present the results of the Volterra-based inverse series for acoustic 
scattering by a finite width well or barrier. Again the reflection and transmission amplitudes can be 



obtained analytically, as can the various terms in the power series for the potential. One can show that 

V 

r{k) = r -=^= . (284) 

(2 - Fo) + 2z Vr^cot^TTT^) 



le 



t{k) = ^ I \ zrr{k) (285) 
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In the case of Vq < 0, one can have any finite value for the magnitude of the interaction (corresponding 
to any finite increase in the velocity of sound in the medium). In the case of a barrier, 0 < Vq< 1 ; 
otherwise one encounters an infinite (Vq=1)ot pure imaginary (Vq>1) velocity of sound. It follows that 

Ki2k)= sin^ Vr^y^ (286) 



The (z) is then 



VAz) = 7-^ I dk ^^-^^ 21g'Ha-2z) (287) 



which is recognized as the Fourier transform of the sinc-fUnction. This is well known to be a square well 
or barrier: 



(288) 



= 0, all other z (289) 



[0129] Rearranging, we find that the region where K is nonzero is 



^min < < -^max (290) 



z 

min 



For a barrier, 0 < < 1 and the first order result has a higher barrier than the true one. For a well, Kq < 
0 and the first order result is shallower than the true one. Thus, although the first order result has the 
correct analytical form of a square well or barrier, it has in correct width and height (or depth). However, 
the explicit form of the result is such that it is trivial to obtain the exact potential firom ^minand z^^a- It can 
be seen that 
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and 



« = ^max + ^min (294) 



2z . Iz 

min max 

a = 



These exact, analytical expressions are found to work very well in computational studies as well Thus, 



it is not necessary to evaluate the beyond y = 1 in order to obtain the exact parameters for a square 

well or barrier interaction. Even so, these higher order terms an also be evaluated analytically. 
[0130] These results can again be compared to those obtained using the Fredhohn-based Bom-Neumann 
inverse scattering series. Razavy [4C] has also obtained an expression for the term. In fact, the result 
is of the form of an infinite series, so a closed expression has not been possible. This also prevented him 
from obtaining higher order corrections. However, the structure manifested at the first order is not a 
simple square well but rather an infinite sequence of steps of decreasing magnitude. Razavy does not 
consider the convergence of the series. Despite these qualitatively incorrect features, it is never-the-less 
possible to use Razavy* s result to determine the square interaction parameters exactly. This is because the 
terms in the infinite series permit one to obtain the correct and a-parameter fi-om the first of the infinite 
series of steps. However, because the Fredhohn-based inversion produces unphysical artifacts that are 
absent firom the Volterra-based results, the latter provides a more robust firamework for an inversion when 
one has an interaction that does not yield an explicit formula for the various terms in the series. 
[0131] It is remarkable that the Volterra-based inverse scattering series for both of these simple potentials 
is able to provide either the exact answer or the exact fimctional form of the interaction with only the first 
order term. Furthermore, the fat that all higher order terms can be evaluated analytically is very usefiil. 
We stress that these results are consequences of the fact that the Volterra-based inversion makes use of 
both the reflection and transmission information. 
XVI. CONCLUSIONS AND FUTURE WORK 

[0132] In this portion of the application, we have used the fat that the acoustic scattering Lippmann- 
Schwinger integral equation (in ID) involving the causal (or anti-causal) Green's fimction can be 
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4 



t 



renonnalized to write it as a Volterra integral equation. Such equations possess the best possible 
convergence behavior under Bom-Neumann iteration. Furthermore, for a wide class of interactions (local, 
differential, or non-local but with the Volterra property), the auxiliary transition operator also possesses 
the Volterra property. Consequently, the inverse acoustic scattering series obtained by reverting the 

Voherra-based series in terms of Vj also converges absolutely and uniformly for all | /| < «>. This does 

not, of course, ensure that the rate of convergence is conveniently rapid. It is well-known that an 
absolutely convergent series can be rearranged or grouped in any manner without affecting its 
convergence [Kaplan [16C]]. Of course, this is not true for divergent or conditionally convergent series. 
In the case of seismic scattering, one may expect the changes in the velocity of sound to be modeled 
reasonably by piece-wise constant interactions sine the distance over which there can be large changes 
should be small compared to distances over which the sonic speed changes less rapidly. 
[0133] Our results show that a Volterra-based inversion can be done as a single comprehensive task, 
provided one has both the reflection and transmission amplitudes as functions of k. Indeed, all ID 
scattering problems that can be formulated in a Lippmann-Schwinger framework have now been shown 
to be invertible, given the r{k) and ^(A:). In subsequent work we shall consider this approach for scattering 
in higher dimensions as well. The implications for various applications such as medial imaging, seismic 
exploration, non-destructive testing, etc, are undercurrent study and results will be reported as they are 
obtained. 

[0134] By appropriate use of Equation (194), we have been able to express the Volterra-based inversion 
in a form that requires only r{k) as input, rather than both r{k) and t(ky This is an important reduction in 
the experimental data required to apply our approach. It has been pointed out to us that for the Dirac delta 
interaction, an approach based on the Heitler damping relation also yields the exact result [17C]. The 
approach is for evaluating a first order approximation only. A complete discussion of the relation to the 
present approach will be given elsewhere. 
APPENDIX C 

[0135] In this Appendix we give a few more details regarding the Fredhohn solution of Equation (28), 

P,(Z) = e"^ + f*^^-^';^(z,z")/7,(z") (296) 

where K{z, z") is defined in Equations (201)-(202). The solution may be written as 
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Pk = ^ + 



D(z,z') 



(297) 



where 



i+oo 
fliz/s:(z,z) + 



2! -to 



J»+oo r+00 
0 Jo 



,K(z\z) K{z\z')j 



(298) 



and 



( 

D(z,z')^ k^^{z,z')-[eyf [J dz^'d^X 



K{z,z') K{z,z'') 
{K(z'\z') K{z'\z'')) 



(299) 



[0136] Note that K{ZyZ) vanishes so long as V{z) is not too singular. Therefore, all diagonal terms in the 
determinants appearing in Equation (298) for 2) vanish. All other terms vanish, as discussed in the text 
above, sine they are of the form TrilC), Consequently, Z) = 1. Now consider the integral 

dz' D{zy z')e'^ . We assume, for simphcity and convenience, that the potential has compact support 

on the domain [0, Z ], and that it is bounded. For any value of k and we conclude that k^fK(Zy z^< |a|, 
where a is some finite number. By Hadamards theorem [6, 14], the value of an nth order determinant 
formed from such elements is bounded by |fl|" ii^^l Then the nth term, say t„, in Equation (299) is bounded 
by 



1 



a 



n 



n 



n/2 



(300) 



[0137] Using Stirling's approximation, one has that 



a 



n 



(301) 



By the root test (Kaplan [16]), we see that 



1. ( \^^^ 
limU„ I = lim 



a 



\ 



= 0 



(302) 
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The r adius o f c onvergence is one d i vided b y t his 1 imit s o w e c onciude t hat t he s eries for D(z,z^ 
converges absolutely independent of the strength of the coupling parameter, ^'or the value of ifc. 
[0138] The most robust treatment of the inverse acoustic scattering problem is that based on the reversion 
of the Bom-Neumann series solution of the Lippmann- Schwinger equation. An important issue for this 
approach to inversion is the radius of convergence of the Bom-Neumann series for Fredholm integral 
kernels, and especially for acoustic scattering for which the interaction depends on the square of the 
frequency. By contrast, it is well known that the Bom-Neumann series for the Volterra integral equations 
in quantum scattering are absolutely convergent, in-dependent of the strength of the coupling 
characterizing the interaction. The transformation of the Lippmann-Schwinger equation from a Fredhohn 
to a Volterra structure by renormalization has been considered previously for quantum scattering 
calculations and electromagnetic scattering. In this portion of the application, we employ the 
renormalization technique to obtain a Volterra equation framework for the inverse acoustic scattering 
series, proving that this series also converges absolutely in the entire complex plane of coupling constant 
and frequency values. The present results are for acoustic scattering in one dimension but the method is 
general. The approach is illustrated by applications to two simple one dimensional models for acoustic 
scattering. 
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3D ACOUSTIC SPHERICAL INTERACTION AND 
NON-SPHERICAL INTERACTION 

[0139J The following equation is the general acoustic case: 



kr 



(303) 



INTRODUCTION - Ricetti Functions 



[01401 Let = krh^ mdJ, = krJi for m = 'A and ^ = 1 for the acoustic case. 



(304) 



(305) 



where Hi = N,^iJ, . Then T, = VP,l (r) 



t, = v^vg;,j, 



so that we can derive the following 



(306) 



(307) 



(308) 



(309) 
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T, = f[\-k\j„){H;,\T] 



(310) 



(311) 



(312) 



P-:(r) = AnY,i%[k)Y,^{r)P:M 

If 



(313) 



where in is data from spectra. 



[0141] At r outside the range of V(r), we have: 

47r 



(314) 



[0142] We now define: 



Data=- 



47r 



Y,{2U\)p{k-r]H:{kr)T^'^ 



(316) 



(317) 



where 



(318) 



(319) 
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and 



(320) 



7j(2) ^ 



(2) 



1+ kTf 



(2) 



(321) 



which yields 



1- 



(2) 



(2) 



(1) 



\+kT, 



(2) 



(322) 



Data= - — Z (2/ + 1) (k ■ r)H; {kr) ' 



I 



(l+kT,^'^) 



(323) 



[0143] Then to first order, F, = T,, so that 



4;r 



Data= - — E (2/ + 1) ^/ • r)H; (kr) 



(1) 



1/ 



(i+^p;p>) 



(324) 



where 



fO) = /j \v\j ) 



(325) 



Ik 



[0144] The higher order corrections are given by: 



/ \ / \ 



(326) 



(327a...) 



The higher order corrections converge absolutely and uniformly allowing analysis of inverse 
scattering spectral data to any level of accuracy. 
Additional Notation for the Non-Spherical Case 
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kr 



- kY, p/-V''^;(^>)./,(ArjK(/;/l/''//''|r')p;(r//''|/'//>') 



(329) 



(330) 



(331) 



diagonal matrix elements are given by: 



T= T 



T=V+ VG T 



(332a-c) 



(333a-b) 



(334a-c) 
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[0145] We need to use the expression GJ^ = Gq^^ + O 



In 



(335) 



If 



(336a-b) 



= V+VGoJ+VOT 
= V[l+OT]+VGoJ 

r = f[i + or] 
r= F + FGo^r 



(337a-c) 



(338) 



i"m" 



(339) 



Solving for r^^^ and plugging into equation (339) for J^'^ (///| /' //' | /r) . Then we take a first order 



treatment and higher order terms again .result firom 



V = -V G V 



(340) 
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where Gq^ is as described previously leading to a data construction shown below: 



(341) 



[0146] Now we can solve the case to the first order where T = V, impUes 



(342) 



(343) 



Now V^^^^[lju\ r 1 k) can be expressed in terms of a set of calculable coefficients with a smaller 

number of unknown and therefore, reduced computational expense for achieving a given degree of 
accuracy in the calculated data verses the actual data. Thus, we can write: 



Kir) = 1 YjV,{lM\r) 



(344) 



[0147] Get Clebsch-Gordan coefficients and a subset of basic integrals, which give rise to a formal 
description of 7'^'^(A^) as follows: 



(1) 



(2) 



'^'\k) = f}^\k) - kr'\k)f^'\k) 



(2) 



(2) 



'^^\k)=[l 



+ kT 



(2) 



(i)]"'f 



(2) 



T^'Kk) = f''\k) 1- ^1+ kT''\k)Y T''\k) 



-1 



(345a-d) 
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We can now define 



V}'Kk) = T^'Hk) (346) 



and we get 



T^'Kk) = V''\k)[u kV''\k)y (347) 



This give us 



T^HWmW = Z K<'>(///|/"//"|^) fl+ kV^'\k)Y\ (348) 



[0148] This expression then is used in the expression for the scatted pressure wave. One determines 
VAr), Then by using the full 3D G^q^ operator in the coordinate representation, we obtain the 

higher order corrections. They will be easier to determine because of their form than to calculate 



or determine Fj (r ) . These equations have utility in ultrasound medical imaging and in sonar for 

submarines and ships to remove reverberations and sonic clutter. 

ELECTROMAGNETIC INVERSE SCATTERING 



vf ) = ^ Z Wfj-) l', (350) 



r/;.fr,^') = ri,.,(r,r')+ 5.1,..,W5i,,„(r') (351) 
[0149] For the case of A = e and = e, we have 

fi...(r,r') = e[u^j ikr')W^;^ {kr) - u'j {h-)w^^;^ {kr')), r'< r 
= 0, r > r 
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[0150] For the case of A, = M and X'" = M, we have 



(353) 



= 0, r'> r 



(354) 



[0151] For the case of X = 0 and X'" = 0, we have 



(355) 



IT ' 

= 0, r'> r 



(356) 



rr 



(357) 



and we have 



= 0, r'> r 



(358) 



= - i-^^-^^/)]" |[.,(^)wy>(A.-)] 



(359) 
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(360) 



where 



[0152] Then we can show that 



and 



Soir) = U{kr) + [ ^0(^0 = S ^oy W 



(361) 



(362) 



(363) 



(364) 



(365) 



(366) 



j=o 



(367) 



Now asymptotically, for /• > r^, the scattering amplitude is obtained from s{r) ■ Thus, one 
determines ;7(]^^J^(,-) from the scattering. Then, one obtains (r) from the spectra to start 



the iterative process. 



SUBTRACTION TECHNIQUE APPROACH TO 
THE INVERSE SCATTERING PROBLEM 
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[0153] Let 



T=V+KT (368) 



where 



K = VG;, (369) 
GL = a. + G.. (370) 



so = JT, + where 



^1 = ^<^o* (371) 



^2 = f^<^o* (372) 
[0154] We now define the following operators as follows: 

r, = i+js:,r, (373) 

r, = (i+jr,)-» (374) 

Then, we can show 

T= r,V+ T,K,T= (1- rao,)'V(l+ G,j) (375) 



So we define the effective potential V as follows: 



V = {l-VOj-'v (376) 

[01551 It can be shown that V is non-local even if Fis local. Therefore, we can show that 

T= V^VG^.T (377) 



We now require that V is local (and also that Vis local). Solving for V , we obtain 

F=r(l+Go,r)"' (378) 
Now we perform a power series expansion of the fraction as follows: 

{uG„r)-' = £(-G„,rr 



«=0 ^ ' 



to obtain the following equation: 
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CO 



n=0 



(380) 



00 



7=1 



(381) 



and T is first order in A, so 



CO 



00 



>t=0 



(382) 



[0156] Then we can write 



(383) 



(384) 



Since F, must be local in coordinate space, it is completely determined by the matrix elements 

(- mk) (Newton, Scattering Theory of Waves and Particles, Springer- Verlag, New York, 1 982, 



chpt. 20). Additionally, off-shell element of are easily gotten as {k'\V^\k) so V can be 



constructed. The physical potential is obtained from equation (376): 



{\-vo,,)v = v 



(385) 



so 



(386) 



[0157] We can evaluate from -k,k matrix elements to solve for F (2 A:) and Fourier invert. 



v=v{\^o^,v) 



-1 



(387) 



The fact that V is local will ensure that V is also local, since Oq* is a separable operator. This is 
an alternative way to obtain the same inversion as before. The only iteration is for the Volterra 
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kernel equation V or T of equation (377). We also note that 



r= F+ FG^F (388) 



The operator Q is Volterra so for V local, T is also Volterra. 
[0158] Now evaluating the equations in momentum space, we have 

{\-vo,,)v = v 



But if we set k=0, then we have 



that is the average of V and Vztq the same. Therefore, we have 



from which V{z) can be obtained by inverse Fourier transformation. 
[0159] Alternatively, 



(389) 



(l + ik7iV\ k){k |)f = V (390) 
V{2k) + ik7rV{2k)V(0) = V{2k) (391) 



ViO) = F(0) (392) 



V{2k) = V(2k)[l - ikTiViO)] (393) 



V(2k) + ik7rV{2k)V{0) = V{2k) (394) 
V(2k) = ^ ^ 



1 - ikT^m (395) 



and this expression is correct since V{0) is essentially t{k). 

MISCELLANEOUS VOLTERRA INVERSE SCATTERING RESULTS 

1 ITT 

= Go,-—\k){k\ (396) 



E-K+ie °* k 



so we now have 
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1= rdr'\r){k' 



(398) 



S(E-K)= I'jk'S^E - k'^)\k'){k' 



(399) 



E'= k'^ 
dE'= 2k' dk' 



(400a-c) 



dk'= 



1 



2^fE' 



dE' 



S{E-K)= l^'^dk'S{E-k'')\k'){k'\+ [jk'S{E-k'^)\km (401) 



[0160] Now let -II = k' and dii = - A', then we can derive 



S(E-K)= [yuS{E - u% u){- u\ (402) 



<5(^ - JiT) = ~l\ ArX^I + I - At I] (403) 



- /;r(5(£ -K)^- ^[\k){k\ + \- k){- k\] (404) 



Therefore, we can write 



P in 

^0* = — (405) 



The principle leading term is part of a Green's function. 

GENERAL APPROACH TO VOLTERRA-BASED 
INVERSION QUANTUM SCATTERING CASE 

[01611 We start with 



T=V^VG',J (406) 
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ft' ^ 



m 



(407) 



We can now separate Gq^ into a Volterra kernel plus a separable Fredholm kernel: 



where 



= ^Ok + ^k 



(408) 



m 



(409) 



where 



{^\y,Jik)=yJr)j,ikr) 



and with hl^{kr) = «/(^) + //,(^) , 



(410) 



(411) 



^0* = - 



2 X X [I X^y/* 



I m 



Y,J ik){Y,jiu,\ , r'<r 



(412) 



= 0, r'> r 

[0162] Using equation (407) and equation (410), we can write 



T=V^ VOJ+ VGqJ 



(413) 



(414) 



We can now separate the Volterra and Fredholm pieces of equation (414) by defining f by 



7= f[l+ Oj] 



(415) 



T=V^-VG^J 



(416) 



To generate V, we expand it in orders oif and require the coefficient of each ^ vanish separately. 
Thus, 
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GO 



(417) 



and 



AT (418) 



[0163] Then by equation (12) 



00 00 



XT=Y, X% + X '^'^y<^o*'^7' (419) 



Therefore 



(420) 



and 



^> = -^y.i<^o.^i. y>l (421) 



[0164] We shall assume that all backscattering elements of T are known, z.e,, all elements of 
(- ^1 J] k) are known. We now define a local operator Fj (r ) such that 

(- ^m^) - ^ (422) 
If (- k\ T\ k) is measured (known), then (f) results fi"om the inverse Fourier transform as shown 



below: 



2ldke-''''{- k\T\k) = ^ p(2^) {^•e^'*(^--V,(r') (423) 



which lead to the following: 



F,(f) = 2ldke-'"'''{- k\T\k) (424) 
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[01651 We note that once (r ) is known, then we can compute its matrix elements in any other 



representation. To determine (f) , we must solve the Fredholm integral equation set forth in 



equation (415). Using equation (410), we have 



v, = v, 



I m 



(425) 



The general matrix element of and V^ , in the angular momentum, radial basis are given below: 
{yiJ,,\vK'Jn) = K j^^'^:WyJ^)^,(^,^K.(^")7/.*(^*) (426) 



But we also have 



J') ^ sir - Py^ir) 



leading to the following 



(427) 



(428) 



[0166] The same is true for the V, -matrix 



(429) 



Then equation (425) yields 



^l"m"Jl"k ^1 l-^'/w'^/'A / - {j^l"m"jl"k ^1 ^im'Jl' 



2mk T-^ 



2 2j ^ {^l"m"Jl"k 



n 



V, 



^ImJlk K^lm^lk 



I m 



(430) 



[0167] This gives a set of linear, inhomogeneous algebraic equations which can be solved for 



Yr„,''ji"k 1^ lUi'k ) • Next, we form 
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2 

= (- k\v,\k) 



(431) 



Then we can show that 



(432) 



[0168] Now to generate (f) , we form the following 









/X/ 















(- 



/ —•I r>w l-w — • \ 

= - jc/r |c/F'e**"T,.,(r)Go,(r,F')^(rOe'*'" 



(433) 



(434) 



and we also form the following 



Vj{r) = 2ldke'^'''V.{2k) 



(435) 



This gives a complete, absolutely convergent scheme, combined with a solution of a separable-kemel 
Fredhohn equation of the second kind (inhomogeneous). All the mathematic proofs and derivations 
are rigorous. 



[0169] One additional step can be to expand and Fj in spherical harmonics as follows: 



(436) 



(437) 



Now we can define 



(438) 



and 
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V,{lrn\rrk)= [ drr^j,.{kr)j,{kr)Vi{im\r) (439) 
V,{lm\rrk)= [drr^J,.{/cr)j,{kr)V,{Im\r) (440) 

V^{lm\l" r k) = f drT^h,.(kr)j,,(kr)vXlm\r) (441) 



From these results, we get the following 



X C(/' m' / w I /• ' m")v^(lm\r r k) = c{V ni lrn\ V 'Tri% (/ w | V'Vk) 



Im Im 



^ SEE C{lmlrn\ /' ' w' (/m | /' ' lk)c[lmV w ' 1 /' m'jFj (/ ' m ' | //' A;) ^"^^^ 



Then we want to solve for the elements , because there are fewer 

elements to calculate than the matrix elements (l^«„..yy„^ \Vy \Yj,^j)*f^ ) • 

[01 70] All r eferences cited h erein a re i ncorporated b y r eference. W hile t his i nvention h as b een 
described fully and completely, it should be understood that, within the scope of the appended claims, 
the invention may be practiced otherwise than as specifically described. Although the invention has 
been disclosed with reference to its preferred embodiments, from reading this description those of skill 
in the art may appreciate changes and modification that may be made which do not depart firom the 
scope and spirit of the invention as described above and claimed hereafter. 
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